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TECHNICAL FIELD 
This invention relates to using molecular and evolutionary techniques to identify 
polynucleotide and polypeptide sequences corresponding to evolved traits that may be 
relevant to human diseases or conditions, such as unique or enhanced human brain 
20 functions, longer human life spans, susceptibility or resistance to development of 
infectious disease (such as AIDS and hepatitis C), susceptibility or resistance to 
development of cancer, and aesthetic traits, such as hair growth, susceptibility or 
resistance to acne, or enhanced muscle mass. 



25 BACKGROUND OF THE INVENTION 

Humans differ from their closest evolutionary relatives, the non-human primates 
such as chimpanzees, in certain physiological and functional traits that relate to areas 
important to human health and well-being. For example, (1) humans have unique or 
enhanced brain function (e.g., cognitive skills, etc.) compared to chimpanzees; (2) 

30 humans have a longer life-span than non-human primates; (3) chimpanzees are resistant 
to certain infectious diseases that afflict humans, such as AIDS and hepatitis C; (4) 
chimpanzees appear to have a lower incidence of certain cancers than humans; (5) 
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chimpanzees do not suffer from acne or alopecia (baldness); (6) chimpanzees have a 
higher percentage of muscle to fat; (7) chimpanzees are more resistant to malaria; (8) 
chimpanzees are less susceptible to Alzheimer's disease; and (9) chimpanzees have a 
lower incidence of atherosclerosis. At the present time, the genes underlying the above 

5 human/chimpanzee differences are not known, nor, more importantly, are the specific 
changes that have evolved in these genes to provide these capabilities. Understanding the 
basis of these differences between humans and our close evolutionary relatives will 
provide useful information for developing effective treatments for related human 
conditions and diseases. 

10 Classic evolution analysis, which compares mainly the anatomic features of 

animals, has revealed dramatic morphological and functional differences between human 
and non-human primates; yet, the human genome is known to share remarkable sequence 
similarities with that of other primates. For example, it is generally concluded that 
human DNA sequence is roughly 98.5% identical to chimpanzee DNA and only slightly 

15 less similar to gorilla DNA. McConkey and Goodman (1997) TIG 13:350-351. Given 
the relatively small percentage of genomic difference between humans and closely related 
primates, it is possible, if not likely, that a relatively small number of changes in genomic 
sequences may be responsible for traits of interest to human health and well-being, such 
as those listed above. Thus, it is desirable and feasible to identify the genes underlying 

20 these traits and to glean information from the evolved changes in the proteins they encode 
to develop treatments that could benefit human health and well-being. Identifying and 
characterizing these sequence changes is crucial in order to benefit from evolutionary 
solutions that have eliminated or minimized diseases or that provide unique or enhanced 
functions. 

25 Recent developments in the human genome project have provided a tremendous 

amount of information on human gene sequences. Furthermore, the structures and 
activities of many human genes and their protein products have been studied either 
directly in human cells in culture or in several animal model systems, such as the 
nematode, fruit fly, zebrafish and mouse. These model systems have great advantages in 
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being relatively simple, easy to manipulate, and having short generation times. Because 
the basic structures and biological activities of many important genes have been 
conserved throughout evolution, homologous genes can be identified in many species by 
comparing macromolecule sequences. Information obtained from lower species on 

5 important gene products and functional domains can be used to help identify the 

homologous genes or functional domains in humans. For example, the homeo domain 
with DNA binding activity first discovered in the fruit fly Drosophila was used to identify 
human homologues that possess similar activities. 

Although comparison of homologous genes or proteins between human and a 

10 lower model organism may provide useful information with respect to evolutionarily 

conserved molecular sequences and functional features, this approach is of limited use in 
identifying genes whose sequences have changed due to natural selection. With the 
advent of the development of sophisticated algorithms and analytical methods, much 
more information can be teased out of DNA sequence changes. The most powerful of 

15 these methods, "K A /Ks" involves pairwise comparisons between aligned protein-coding 
nucleotide sequences of the ratios of 

nonsvnonvmous nucleotide substitutions per nonsvnonvmous site (Ka) 
synonymous substitutions per synonymous site (Ks) 

20 

(where nonsynonymous means substitutions that change the encoded amino acid and 
synonymous means substitutions that do not change the encoded amino acid). "K A /Ks- 
type methods" includes this and similar methods. These methods have been used to 
demonstrate the occurrence of Darwinian molecular-level positive selection, resulting in 

25 amino acid differences in homologous proteins. Several groups have used such methods 
to document that a particular protein has evolved more rapidly than the neutral 
substitution rate, and thus supports the existence of Darwinian molecular-level positive 
selection. For example, McDonald and Kreitman (1991) Nature 351 :652-654 propose a 
statistical test of neutral protein evolution hypothesis based on comparison of the number 

30 of amino acid replacement substitutions to synonymous substitutions in the coding region 
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of a locus. When they apply this test to the Adh locus of three Drosophila species, they 
conclude that it shows instead that the locus has undergone adaptive fixation of 
selectively advantageous mutations and that selective fixation of adaptive mutations may 
be a viable alternative to the clocklike accumulation of neutral mutations as an 

5 explanation for most protein evolution. Jenkins et al. (1995) Proa R. Soc. Lond. B 
261:203-207 use the McDonald & Kxeitman test to investigate whether adaptive 
evolution is occurring in sequences controlling transcription (non-coding sequences). 

Nakashima et al. (1995) Proc. Natl Acad Sci USA 92:5606-5609, use the method 
of Miyata and Yasunaga to perform pairwise comparisons of the nucleotide sequences of 

10 ten PLA2 isozyme genes from two snake species; this method involves comparing the 
number of nucleotide substitutions per site for the noncoding regions including introns 
(K N ) and the K A and K s . They conclude that the protein coding regions have been 
evolving at much higher rates than the noncoding regions including introns. The highly 
accelerated substitution rate is responsible for Darwinian molecular-level evolution of 

1 5 PLA2 isozyme genes to produce new physiological activities that must have provided 
strong selective advantage for catching prey or for defense against predators. Endo et al 
(1996) Mol Biol Evol 13(5):685-690 use the method of Nei and Gojobori, wherein d N is 
the number of nonsynonymous substitutions and d s is the number of synonymous 
substitutions, for the purpose of identifying candidate genes on which positive selection 

20 operates. Metz and Palumbi (1996) Mol Biol Evol 13(2):397-406 use the McDonald & 
Kreitman test as well as a method attributed to Nei and Gojobori, Nei and Jin, and 
Kumar, Tamura, and Nei; examining the average proportions of P n , the replacement 
substitutions per replacement site, and P s , the silent substitutions per silent site, to look 
for evidence of positive selection on bindin genes in sea urchins to investigate whether 

25 they have rapidly evolved as a prelude to species formation. Goodwin et al (1996) Mol 
Biol Evol 13(2):346-358 uses similar methods to examine the evolution of a particular 
murine gene family and conclude that the methods provide important fundamental 
insights into how selection drives genetic divergence in an experimentally manipulatable 
system. Edwards et al (1995) use degenerate primers to pull out MHC loci from various 
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species of birds and an alligator species, which are then analyzed by the Nei and Gojobori 
methods (d N : ds ratios) to extend MHC studies to nonmammalian vertebrates. Whitfield 
et aL (1993) Nature 364:713-715 use Ka/Ks analysis to look for directional selection in 
the regions flanking a conserved region in the 57? F gene (that determines male sex). They 
5 suggest that the rapid evolution of SR Y could be a significant cause of reproductive 
isolation, leading to new species. Wettsetin et aL (1996) MoL BioL EvoL 13(l):56-66 
apply the MEGA program of Kumar, Tamura and Nei and phylogenetic analysis to 
investigate the diversification of MHC class I genes in squirrels and related rodents. 
Parham and Ohta (1996) Science 212:61-1 A state that a population biology approach, 
10 including tests for selection as well as for gene conversion and neutral drift are required 
to analyze the generation and maintenance of human MHC class I polymorphism. 
Hughes (1997) MoL Biol. EvoL 14(1): 1-5 compared over one hundred orthologous 
immunoglobulin C2 domains between human and rodent, using the method of Nei and 
Gojobori (d N : d s ratios) to test the hypothesis that proteins expressed in cells of the 
15 vertebrate immune system evolve unusually rapidly. Swanson and Vacquier (1998) 
Science 281:710-712 use d N : d s ratios to demonstrate concerted evolution between the 
lysin and the egg receptor for lysin and discuss the role of such concerted evolution in 
forming new species (speciation). 

Due to the distant evolutionary relationships between humans and these lower 
20 animals, the adaptively valuable genetic changes fixed by natural selection are often 
masked by the accumulation of neutral, random mutations over time. Moreover, some 
proteins evolve in an episodic manner; such episodic changes could be masked, leading to 
inconclusive results, if the two genomes compared are not close enough. Messier and 
Stewart (1997) Nature 385:151-154. In fact, studies have shown that the occurrence of 
25 adaptive selection in protein evolution is often underestimated when predominantly 
distantly related sequences are compared. Endo et al. (1996) MoL BioL EvoL 37:441- 
456; Messier and Stewart (1997) Nature 385:151-154. 

Molecular evolution studies within the primate family have been reported, but 
these mainly focus on the comparison of a small number of known individual genes and 
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gene products to assess the rates and patterns of molecular changes and to explore the 
evolutionary mechanisms responsible for such changes. See generally, Li, Molecular 
Evolution, Sinauer Associates, Sunderland, MA, 1997, Furthermore, sequence 
comparison data are used for phylogenetic analysis, wherein the evolution history of 
5 primates is reconstructed based on the relative extent of sequence similarities among 
examined molecules from different primates. For example, the DNA and amino acid 
sequence data for the enzyme lysozyme from different primates were used to study 
protein evolution in primates and the occurrence of adaptive selection within specific 
lineages. Malcolm et al (1990) Nature 345:86-89; Messier and Stewart (1997). Other 

10 genes that have been subject to molecular evolution studies in primates include 

hemoglobin, cytochrome c oxidase, and major histocompatibility complex {MHC). Nei 
and Hughes in: Evolution at the Molecular Level, Sinauer Associates, Sunderland, MA 
222-247, 1991; Lienert and Parham (1996) Immunol Cell Biol 74:349-356; Wu et al 
(1997) J. Mol Evol 44:477 -491. Many non-coding sequences have also been used in 

15 molecular phylogenetic analysis of primates. Li, Molecular Evolution, Sinauer 

Associates, Sunderland, MA 1997. For example, the genetic distances among primate 
lineages were estimated from orthologous non-coding nucleotide sequences of beta-type 
globin loci and their flanking regions, and the evolution tree constructed for the 
nucleotide sequence orthologues depicted a branching pattern that is largely congruent 

20 with the picture from phylogenetic analyses of morphological characters. Goodman et al. 
(1990)/. Mol Evol 30:260-266. 

Zhou and Li (1996) Mol Biol Evol 13(6):780-783 applied K A /K S analysis to 
primate genes. It had previously been reported that gene conversion events likely have 
occurred in introns 2 and 4 between the red and green retinal pigment genes during 

25 human evolution. However, intron 4 sequences of the red and green retinal pigment 
genes from one European human were completely identical, suggesting a recent gene 
conversion event. In order to determine if the gene conversion event occurred in that 
individual, or a common ancestor of Europeans, or an even earlier hominid ancestor, the 
authors sequenced intron 4 of the red and green pigment gene from a male Asian human, 
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a male chimpanzee, and a male baboon, and applied K A /Ks analysis. They observed that 
the divergence between the two genes is significantly lower in intron 4 than in 
surrounding exons, suggesting that strong natural selection has acted against sequence 
homogenization. 

5 Wolinsky et ai (1996) Science 272:537-542 used comparisons of nonsynonymous 

to synonymous base substitutions to demonstrate that the HIV virus itself (i.e., not the 
host species) is subject to adaptive evolution within individual human patients. Their 
goal was simply to document the occurrence of positive selection in a short time frame 
(that of a human patient's course of disease). Niewiesk and Bangham (1996) JMol Evol 
10 42:452-458 used the D n /D s approach to ask a related question about the HTLV-1 virus, 
i.e., what are the selective forces acting on the virus itself. Perhaps because of an 
insufficient sample size, they were unable to resolve the nature of the selective forces. In 
both of these cases, although K A /K s -type methods were used in relation to a human virus, 
no attempt was made to use these methods for therapeutic goals (as in the present 
1 5 application), but rather to pursue narrow academic goals. 

As can be seen from the papers cited above, analytical methods of molecular evolution to 
identify rapidly evolving genes (K A /K s -type methods) can be applied to achieve many 
different purposes, most commonly to confirm the existence of Darwinian molecular- 
level positive selection, but also to assess the frequency of Darwinian molecular-level 
20 positive selection, to understand phylogenetic relationships, to elucidate mechanisms by 
which new species are formed, or to establish single or multiple origin for specific gene 
polymorphisms. What is clear is from the papers cited above and others in the literature 
is that none of the authors applied K A /K s -type methods to identify evolutionary solutions, 
specific evolved changes, that could be mimicked or used in the development of 
25 treatments to prevent or cure human conditions or diseases or to modulate unique or 

enhanced human functions. They have not used K A /K S type analysis as a systematic tool 
for identifying human or non-human primate genes that contain evolutionarily significant 
sequence changes and exploiting such genes and the identified changes in the 
development of treatments for human conditions or diseases. 
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The identification of human genes that have evolved to confer unique or enhanced 
human functions compared to homologous chimpanzee genes could be applied to 
developing agents to modulate these unique human functions or to restore function when 
the gene is defective. The identification of the underlying chimpanzee (or other non- 

5 human primate) genes and the specific nucleotide changes that have evolved, and the 
further characterization of the physical and biochemical changes in the proteins encoded 
by these evolved genes, could provide valuable information, for example, on what 
determines susceptibility and resistance to infectious diseases, such as AIDS, what 
determines susceptibility or resistance to the development of certain cancers, what 

10 determines susceptibility or resistance to acne, how hair growth can be controlled, and 
how to control the formation of muscle versus fat. This valuable information could be 
applied to developing agents that cause the human proteins to behave more like their 
chimpanzee homologues. 

All references cited herein are hereby incorporated by reference in their entirety. 

15 

SUMMARY OF THE INVENTION 
The present invention provides methods for identifying polynucleotide and 
polypeptide sequences having evolutionarily significant changes which are associated 
with physiological conditions, including medical conditions. The invention applies 
20 comparative primate genomics to identify specific gene changes which may be associated 
with, and thus responsible for, physiological conditions, such as medically or 
commercially relevant evolved traits, and using the information obtained from these 
evolved genes to develop human treatments. The non-human primate sequences 
employed in the methods described herein may be any non-human primate, and is 
25 preferably a member of the hominoid group, more preferably a chimpanzee, bonobo, 
gorilla and/or orangutan, and most preferably a chimpanzee. 

In one preferred embodiment, a non-human primate polynucleotide or polypeptide 
has undergone natural selection that resulted in a positive evolutionarily significant 
change (i.e., the non-human primate polynucleotide or polypeptide has a positive attribute 



not present in humans). In this embodiment the positively selected polynucleotide or 
polypeptide may be associated with susceptibility or resistance to certain diseases or with 
other commercially relevant traits. Examples of this embodiment include, but are not 
limited to, polynucleotides and polypeptides that are positively selected in non-human 

5 primates, preferably chimpanzees, that may be associated with susceptibility or resistance 
to infectious diseases and cancer. An example of a commercially relevant trait may 
include aesthetic traits such as hair growth, muscle mass, susceptibility or resistance to 
acne. An example of this embodiment includes polynucleotides and polypeptides 
associated with the susceptibility or resistance to HIV dissemination, propagation and/or 

10 development of AIDS. The present invention can thus be useful in gaining insight into 
the molecular mechanisms that underlie resistance to HIV dissemination, propagation 
and/or development of AIDS, providing information that can also be useful in discovering 
and/or designing agents such as drugs that prevent and/or delay development of AIDS. 
Specific genes that have been positively selected in chimpanzees that may relate to AIDS 

15 or other infectious diseases are ICAM-1, ICAM-2, ICAM-3, MIP-1-a, CD59 and DC- 
SIGN. 17-p-hydroxysteroid dehydrogenase Type IV is a specific gene has been positively 
selected in chimpanzees that may relate to cancer. 

In another preferred embodiment, a human polynucleotide or polypeptide has 
undergone natural selection that resulted in a positive evolutionarily significant change 

20 (i.e., the human polynucleotide or polypeptide has a positive attribute not present in non- 
human primates). One example of this embodiment is that the polynucleotide or 
polypeptide may be associated with unique or enhanced functional capabilities of the 
human brain compared to non-human primates. Another is the longer life-span of 
humans compared to non-human primates. A third is a commercially important aesthetic 

25 trait (e.g., normal or enhanced breast development). The present invention can thus be 
useful in gaining insight into the molecular mechanisms that underlie unique or enhanced 
human functions or physiological traits, providing information which can also be useful 
in designing agents such as drugs that modulate such unique or enhanced human 
functions or traits, and in designing treatment of diseases or conditions related to humans. 
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As an example, the present invention can thus be useful in gaining insight into the 
molecular mechanisms that underlie human cognitive function, providing information 
which can also be useful in designing agents such as drugs that enhance human brain 
function, and in designing treatment of diseases related to the human brain. A specific 
5 example of a human gene that has positive evolutionarily significant changes when 
compared to non-human primates is a tyrosine kinase gene, the KIAA0641 or 
NM_004920 gene. 

Accordingly, in one aspect, the invention provides methods for identifying a 
polynucleotide sequence encoding a polypeptide, wherein said polypeptide may be 

10 associated with a physiological condition (such as a medically or commercially relevant 
positive evolutionarily significant change). The positive evolutionarily significant change 
can be found in humans or in non-human primates. In a preferred embodiment the 
invention provides a method for identifying a human AATYK polynucleotide sequence 
encoding a human AATYK polypeptide associated with an evolutionarily significant 

15 change. 

For any embodiment of this invention, the physiological condition may be any 
physiological condition, including those listed herein, such as, for example, disease 
(including susceptibility or resistance to disease) such as cancer, infectious disease 
(including viral diseases such as AIDS); life span; brain function, including cognitive 

20 function or developmental sculpting; and aesthetic or cosmetic qualities, such as 
enhanced breast development. 

In one aspect of the invention, methods are provided for identifying a 
polynucleotide sequence encoding a human polypeptide, wherein said polypeptide may be 
associated with a physiological condition that is present in human(s), comprising the 

25 steps of: a) comparing human protein-coding polynucleotide sequences to protein-coding 
polynucleotide sequences of a non-human primate, wherein the non-human primate does 
not have the physiological condition (or has it to a lesser degree); and b) selecting a 
human polynucleotide sequence that contains a nucleotide change as compared to 
corresponding sequence of the non-human primate, wherein said change is evolutionarily 
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significant. In some embodiments, the human protein coding sequence (and/or the 
polypeptide encoded therein) may be associated with development and/or maintenance of 
a physiological condition or trait or a biological function. In some embodiments, the 
physiological condition or biological function may be life span, brain or cognitive 
5 function, or breast development (including adipose, gland and duct development). 

Methods used to assess the nucleotide change, and the nature(s) of the nucleotide change, 
are described herein, and apply to any and all embodiments. In a preferred embodiment, 
the method is a method for identifying a human AATYK polynucleotide sequence 
encoding a human AATYK polypeptide. 
10 In other embodiments, methods are provided that comprise the steps of: (a) 

comparing human protein-coding nucleotide sequences to protein-coding nucleotide 
sequences of a non-human primate, preferably a chimpanzee, that is resistant to a 
particular medically relevant disease state, wherein the human protein coding sequence is 
associated with development of the disease; and (b) selecting a non-human 
15 polynucleotide sequence that contains at least one nucleotide change as compared to the 
corresponding sequence of the human, wherein the change is evolutionarily significant. 
The sequences identified by these methods may be further characterized and/or analyzed 
to confirm that they are associated with the development of the disease state or condition. 
The most preferred disease states that are applicable to these methods are cancer and 
20 infectious diseases, including AIDS, hepatitis C and leprosy. 

In another aspect, methods are provided for identifying an evolutionarily 
significant change in a human brain polypeptide-coding polynucleotide sequence, 
comprising the steps of a) comparing human brain polypeptide-coding polynucleotide 
sequences to corresponding sequences of a non-human primate; and b) selecting a human 
25 polynucleotide sequence that contains a nucleotide change as compared to corresponding 
sequence of the non-human primate, wherein said change is evolutionarily significant. In 
some embodiments, the human brain polypeptide coding nucleotide sequences 
correspond to human brain cDNAs. In preferred embodiments, the human brain 
polypeptide-coding polynucleotide sequence is an AATYK sequence. 



Another aspect of the invention includes methods for identifying a positively 
selected human evolutionarily significant change. These methods comprise the steps of: 
(a) comparing human polypeptide-coding nucleotide sequences to polypeptide-coding 
nucleotide sequences of a non-human primate; and (b) selecting a human polynucleotide 
sequence that contains at least one (i.e., one or more) nucleotide change as compared to 
corresponding sequence of the non-human primate, wherein said change is evolutionarily 
significant. The sequences identified by this method maybe further characterized and/or 
analyzed for their possible association with biologically or medically relevant functions or 
traits unique or enhanced in humans. In preferred embodiments, the human polypeptide- 
coding nucleotide sequence is an AATYK sequence. 

Another embodiment of the present invention is a method for large scale sequence 
comparison between human polypeptide-coding polynucleotide sequences and the 
polypeptide-coding polynucleotide sequences from a non-human primate, e.g., 
chimpanzee, comprising: (a) aligning the human polynucleotide sequences with 
corresponding polynucleotide sequences from non-human primate according to sequence 
homology; and (b) identifying any nucleotide changes within the human sequences as 
compared to the homologous sequences from the non-human primate, wherein the 
changes are evolutionarily significant. In some embodiments, the protein coding 
sequences are from brain. 

In some embodiments, a nucleotide change identified by any of the methods 
described herein is a non-synonymous substitution. In some embodiments, the 
evolutionary significance of the nucleotide change is determined according to the non- 
synonymous substitution rate (K A ) of the nucleotide sequence. Li some embodiments, the 
evolutionarily significant changes are assessed by determining the K A /K S ratio between 
the human gene and the homologous gene from non-human primate (such as 
chimpanzee), and preferably that ratio is at least about 0.75, more preferably greater than 
about 1 (unity) (i.e., at least about 1), more preferably at least about 1.25, more preferably 
at least about 1.50, and more preferably at least about 2.00. In other embodiments, once a 
positively selected gene has been identified between human and a non-human primate 
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(such as chimpanzee or gorilla), further comparisons are performed with other non- 
human primates to confirm whether the human or the non-human primate (such as 
chimpanzee or gorilla) gene has undergone positive selection. 

In another aspect, the invention provides methods for correlating an evolutionarily 

5 significant human nucleotide change to a physiological condition in a human (or 

humans), which comprise analyzing a functional effect (which includes determining the 
presence of a functional effect), if any, of (the presence or absence of) a polynucleotide 
sequence identified by any of the methods described herein, wherein presence of a 
functional effect indicates a correlation between the evolutionarily significant nucleotide 

10 change and the physiological condition. Alternatively, in these methods, a functional 
effect (if any) may be assessed using a polypeptide sequence (or a portion of the 
polypeptide sequence) encoded by a nucleotide sequence identified by any of the methods 
described herein. In a preferred embodiment, the polynucleotide sequence or polypeptide 
sequence is an AATYK polynucleotide sequence or polypeptide sequence. In a more 

15 preferred embodiment, the AATYK polynucleotide sequences are nucleotides 2180-2329 
of SEQ ID NO:14, nucleotides 2978-3478 of SEQ ID NO:14, and nucleotides 3380-3988 
ofSEQIDNO:14. 

The present invention also provides comparison of the identified polypeptides by 
physical and biochemical methods widely used in the art to determine the structural or 

20 biochemical consequences of the evolutionarily significant changes. Physical methods 
are meant to include methods that are used to examine structural changes to proteins 
encoded by genes found to have undergone adaptive evolution. Side-by-side comparison 
of the three-dimensional structures of a protein (either human or non-human primate) and 
the evolved homologous protein (either non-human primate or human, respectively) will 

25 provide valuable information for developing treatments for related human conditions and 
diseases. For example, using the methods of the present invention, the chimpanzee 
ICAM-1 gene was identified as having positive evolutionary changes compared to human 
ICAM-1 . In a three-dimensional model of two functional domains of the human ICAM-1 
protein it can be seen that five of the six amino acids that have been changed in 
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chimpanzees are immediately adjacent to (i.e., physically touching) amino acid residues 
known to be crucial for binding to the ICAM-1 counter-receptor, LFA-1 ; in each case, the 
human amino acid has been replaced by a larger amino acid in the chimpanzee ICAM-L 
Such information allows insight into designing appropriate therapeutic intervention(s). 

5 Accordingly, in another aspect, the invention provides methods for identifying a target 
site (which includes one or more target sites) which may be suitable for therapeutic 
intervention, comprising comparing a human polypeptide (or a portion of the polypeptide) 
encoded in a sequence identified by any of the methods described herein, with a 
corresponding non-human polypeptide (or a portion of the polypeptide), wherein a 

10 location of a molecular difference, if any, indicates a target site. 

In another aspect, the invention provides methods for identifying a target site 
(which includes one or more target sites) which maybe suitable for therapeutic 
intervention, comprising comparing a human polypeptide (or a portion of the polypeptide) 
encoded in a sequence identified by any of the methods described herein, with a 

1 5 corresponding non-human primate polypeptide (or a portion of the polypeptide), wherein 
a location of a molecular difference, such as an amino acid difference, if any, indicates a 
target site. Target sites can also be nonsynonymous nucleotide changes in a positively 
selected polynucleotide identified by any of the methods described herein. In preferred 
embodiments, the target site is a site on a human AATYK polypeptide. 

20 Biochemical methods are meant to include methods that are used to examine 

functional differences, such as binding specificity, binding strength, or optimal binding 
conditions, for a protein encoded by a gene that has undergone adaptive evolution. Side- 
by-side comparison of biochemical characteristics of a protein (either human or non- 
human primate) and the evolved homologous protein (either non-human primate or 

25 human, respectively) will reveal valuable information for developing treatments for 
related human conditions and diseases. 

In another aspect, the invention provides methods of identifying an agent which 
may modulate a physiological condition, said method comprising contacting an agent 
(i.e., at least one agent to be tested) with a cell that has been transfected with a 
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polynucleotide sequence identified by any of the methods described herein, wherein an 
agent is identified by its ability to modulate function of the polynucleotide sequence. In 
other embodiments, the invention provides methods of identifying an agent which may 
modulate a physiological condition, said method comprising contacting an agent (i.e., at 
least one agent) to be tested with a polypeptide (or a fragment of a polypeptide and/or a 
composition comprising a polypeptide or fragment of a polypeptide) encoded in or within 
a polynucleotide identified by any of the methods described herein, wherein an agent is 
identified by its ability to modulate function of the polypeptide. In preferred 
embodiments of these methods the polynucleotide sequence is an evolutionarily 
significant human AATYK polynucleotide sequence. In more preferred embodiments, 
the polynucleotide sequence is nucleotides 2180-2329 of SEQ ID NO:14, nucleotides 
2978-3478 of SEQ ID NO: 14, and nucleotides 3380-3988 of SEQ ED NO: 14. The 
invention also provides agents which are identified using the screening methods described 
herein. 

In another aspect, the invention provides methods of screening agents which may 
modulate the activity of the human polynucleotide or polypeptide to either modulate a 
unique or enhanced human function or trait or to mimic the non-human primate trait of 
interest, such as susceptibility or resistance to development of a disease, such as AIDS. 
These methods comprise contacting a cell which has been transfected with a 
polynucleotide sequence with an agent to be tested, and identifying agents based on their 
ability to modulate function of the polynucleotide or contacting a polyeptide preparation 
with an agent to be tested and identifying agents based upon their ability to modulate 
function of the polypeptide. In preferred embodiments, the polynucleotide sequence is an 
evolutionarily significant human AATYK polynucleotide sequence. In more preferred 
embodiments, the polynucleotide sequence is nucleotides 2180-2329 of SEQ ID NO:14, 
nucleotides 2978-3478 of SEQ ID NO: 14, and nucleotides 3380-3988 of SEQ ID NO:14. 



BRIEF DESCRIPTION OF THE DRAWINGS 
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Figure 1 depicts a phylogenetic tree for primates within the hominoid group. The 
branching orders are based on well-supported mitochondrial DNA phylogenies. Messier 
and Stewart (1997) Nature 385:151-154. 

Figure 2 (SEQ ID NOS:l-3) is a nucleotide sequence alignment between human 
5 and chimpanzee ICAM-1 sequences (GenBank® accession numbers X06990 and 

X86848, respectively). The amino acid translation of the chimpanzee sequence is shown 
below the alignment. 

Figure 3 shows the nucleotide sequence of gorilla ICAM-1 (SEQ ID NO:4). 
Figure 4 shows the nucleotide sequence of orangutan ICAM-1 (SEQ ID NO: 5). 
10 Figures 5(A)-(E) show the polypeptide sequence alignment of ICAM-1 from 

several primate species (SEQ ID NO:6). 

Figures 6(A)-(B) show the polypeptide sequence alignment of ICAM-2 from 
several primate species (SEQ ID NO:7). 

Figures 7(A)-(D) show the polypeptide sequence alignment of ICAM-3 from 
15 several primate species (SEQ ID NO:8). 

Figure 8 depicts a schematic representation of a procedure for comparing 
human/primate brain polynucleotides, selecting sequences with evolutionarily significant 
changes, and further characterizing the selected sequences. The diagram of Figure 8 
illustrates a preferred embodiment of the invention and together with the description 
20 serves to explain the principles of the invention, along with elaboration and optional 

additional steps. It is understood that any human/primate polynucleotide sequence can be 
compared by a similar procedure and that the procedure is not limited to brain 
polynucleotides. 

Figure 9 illustrates the known phylogenetic tree for the species compared in 
25 Example 1 4, with values of b N and b s mapped upon appropriate branches. Values of b N 
and b s were calculated by the method described in Zhang et al (1998) Proc. Natl Acad. 
Set USA 95:3708-3713. Values are shown above the branches; all values are shown 
100X, for reasons of clarity. Statistical significance was calculated as for comparisons in 
Table 5 (Example 14), and levels of statistical significance are as shown as in Table 5. 
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Note that only the branch leading from the human/chimpanzee common ancestor to 
modern humans shows a statistically significant value for b s - b s . 

Figure 10 illustrates a space-filling model of human CD59 with the duplicated 
GPI link (Asn) indicated by the darkest shading. This GPI link is duplicated in 
5 chimpanzees so that chimp CD59 contains 3 GPI links. The three areas of intermediate 
shading in Figure 10 are other residues which differ between chimp and human. 

Figure 1 1 shows the coding sequence of human DC-SIGN (Genbank Acc. No. 
M98457) (SEQ. ID. NO. 9). 

Figure 12 shows the coding sequence of chimpanzee DC-SIGN (SEQ. ID. NO. 

10 10). 

Figure 13 shows the coding sequence of gorilla DC-SIGN (SEQ. ID. NO. 1 1). 
Figure 14A shows the nucleotide sequence of the human AATYK gene. Start and 
stop codons are underlined (SEQ ID NO: 14). 

Figure 14B shows an 1207 amino acid sequence of the human AATYK gene 
15 (SEQEDNO:16). 

Figure 15 A shows an 1806 base-pair region of the chimp AATYK gene (SEQ ID 

NO:17). 

Figure 15B shows an 1785 base-pair region of the gorilla AATYK gene (SEQ ID 
NO: 18). 

20 

DETAILED DESCRIPTION OF THE INVENTION 
The present invention applies comparative genomics to identify specific gene 
changes which are associated with, and thus may contribute to or be responsible for, 
physiological conditions, such as medically or commercially relevant evolved traits. The 
25 invention comprises a comparative genomics approach to identify specific gene changes 
responsible for differences in functions and diseases distinguishing humans from other 
non-humans, particularly primates, and most preferably chimpanzees, including the two 
known species, common chimpanzees and bonobos (pygmy chimpanzees). For example, 
chimpanzees and humans are 98.5% identical at the DNA sequence level and the present 
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invention can identify the adaptive molecular changes underlying differences between the 
species in a number of areas, including unique or enhanced human cognitive abilities or 
physiological traits and chimpanzee resistance to AIDS and certain cancers. Unlike 
traditional genomics, which merely identifies genes, the present invention provides exact 

5 information on evolutionary solutions that eliminate disease or provide unique or 

enhanced functions or traits. The present invention identifies genes that have evolved to 
confer an evolutionary advantage and the specific evolved changes. 

The present invention results from the observation that human protein-coding 
polynucleotides may contain sequence changes that are found in humans but not in other 

1 0 evolutionarily closely related species such as non-human primates, as a result of adaptive 
selection during evolution. 

The present invention further results from the observation that the genetic 
information of non-human primates may contain changes that are found in a particular 
non-human primate but not in humans, as a result of adaptive selection during evolution. 

15 In this embodiment, a non-human primate polynucleotide or polypeptide has undergone 
natural selection that resulted in a positive evolutionarily significant change (i.e., the non- 
human primate polynucleotide or polypeptide has a positive attribute not present in 
humans). In this embodiment the positively selected polynucleotide or polypeptide may 
be associated with susceptibility or resistance to certain diseases or other commercially 

20 relevant traits. Medically relevant examples of this embodiment include, but are not 
limited to, polynucleotides and polypeptides that are positively selected in non-human 
primates, preferably chimpanzees, that may be associated with susceptibility or resistance 
to infectious diseases and cancer. An example of this embodiment includes 
polynucleotides and polypeptides associated with the susceptibility or resistance to 

25 progression from HIV infection to development of AIDS. The present invention can thus 
be useful in gaining insight into the molecular mechanisms that underlie resistance to 
progression from HIV infection to development of AIDS, providing information that can 
also be useful in discovering and/or designing agents such as drugs that prevent and/or 
delay development of AIDS. Commercially relevant examples include, but are not 
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limited to, polynucleotides and polypeptides that are positively selected in non-human 
primates that may be associated with aesthetic traits, such as hair growth, acne or muscle 
mass. 

Positively selected human evolutionarily significant changes in polynucleotide 

5 and polypeptide sequences may be attributed to human capabilities that provide humans 
with competitive advantages, particularly when compared to the closest evolutionary 
relative, chimpanzee, such as unique or enhanced human brain functions. The present 
invention identifies human genes that evolved to provide unique or enhanced human 
cognitive abilities and the actual protein changes that confer functional differences will be 

10 quite useful in therapeutic approaches to treat cognitive deficiencies as well as cognitive 
enhancement for the general population. 

Other positively selected human evolutionarily significant changes include those 
sequences that may be attributed to human physiological traits or conditions that are 
enhanced or unique relative to close evolutionary relatives, such as the chimpanzee, 

15 including enhanced breast development. The present invention provides a method of 
determining whether a polynucleotide sequence in humans that may be associated with 
enhanced breast development has undergone an evolutionarily significant change relative 
to a corresponding polynucleotide sequence in a closely related non-human primate. The 
identification of evolutionarily significant changes in the human polynucleotide that is 

20 involved in the development of unique or enhanced human physiological traits is 

important in the development of agents or drugs that can modulate the activity or function 
of the human polynucleotide or its encoded polypeptide. 

The practice of the present invention employs, unless otherwise indicated, 
conventional techniques of molecular biology, genetics and molecular evolution, which 

25 are within the skill of the art. Such techniques are explained fully in the literature, such 
as: "Molecular Cloning: A Laboratory Manual", second edition (Sambrook et a/., 1989); 
"Oligonucleotide Synthesis" (MJ. Gait, ed., 1984); "Current Protocols in Molecular 
Biology" (F.M. Ausubel et aU eds., 1987); "PCR: The Polymerase Chain Reaction", 
(Mullis et aL, eds., 1994); "Molecular Evolution", (Li, 1997). 
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Definitions 

As used herein, a "polynucleotide" refers to a polymeric form of nucleotides of 
any length, either ribonucleotides or deoxyribonucleotides, or analogs thereof. This term 
5 refers to the primary structure of the molecule, and thus includes double- and single- 
stranded DNA, as well as double- and single-stranded RNA. It also includes modified 
polynucleotides such as methylated and/or capped polynucleotides. The terms 
"polynucleotide" and "nucleotide sequence" are used interchangeably. 

As used herein, a "gene" refers to a polynucleotide or portion of a polynucleotide 
10 comprising a sequence that encodes a protein. It is well understood in the art that a gene 
also comprises non-coding sequences, such as 5' and 3' flanking sequences (such as 
promoters, enhancers, repressors, and other regulatory sequences) as well as introns. 

The terms "polypeptide," "peptide," and "protein" are used interchangeably herein 
to refer to polymers of amino acids of any length. These terms also include proteins that 
1 5 are post-translationally modified through reactions that include glycosylation, acetylation 
and phosphorylation. 

A "physiological condition" is a term well-understood in the art and means any 
condition or state that can be measured and/or observed. A "physiological condition" 
includes, but is not limited to, a physical condition, such as degree of body fat, alopecia 
20 (baldness), acne or enhanced breast development; life-expectancy; disease states (which 
include susceptibility and/or resistance to diseases), such as cancer or infectious diseases. 
Examples of physiological conditions are provided below (see, e.g., definitions of 
"human medically relevant medical condition", "human commercially relevant 
condition", "medically relevant evolved trait", and "commercially relevant evolved trait") 
25 and throughout the specification, and it is understood that these terms and examples refer 
to a physiological condition. A physiological condition may be, but is not necessarily, the 
result of multiple factors, any of which in turn may be considered a physiological 
condition. A physiological condition which is "present" in a human or non-human 
primate occurs within a given population, and includes those physiological conditions 
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which are unique and/or enhanced in a given population when compared to another 
population. 

The terms "human medically relevant condition" or "human commercially 
relevant condition" are used herein to refer to human conditions for which medical or 
5 non-medical intervention is desired. 

The term "medically relevant evolved trait" is used herein to refer to traits that 
have evolved in humans or non-human primates whose analysis could provide 
information (e.g., physical or biochemical data) relevant to the development of a human 
medical treatment. 

1 0 The term "commercially relevant evolved trait" is used herein to refer to traits that 

have evolved in humans or non-human primates whose analysis could provide 
information (e.g., physical or biochemical data) relevant to the development of a medical 
or non-medical product or treatment for human use. 

The term "K A /K s -type methods" means methods that evaluate differences, 

15 frequently (but not always) shown as a ratio, between the number of nonsynonymous 
substitutions and synonymous substitutions in homologous genes (including the more 
rigorous methods that determine non-synonymous and synonymous sites). These 
methods are designated using several systems of nomenclature, including but not limited 

to K A /K S , d N /d s , D N /D S . 

20 The terms "evolutionarily significant change" or "adaptive evolutionary change" 

refers to one or more nucleotide or peptide sequence change(s) between two species that 
may be attributed to a positive selective pressure. One method for determining the 
presence of an evolutionarily significant change is to apply a K A /Ks-type analytical 
method, such as to measure a K A /K S ratio. Typically, a K A /K S ratio at least about 0.75, 

25 more preferably at least about 1 .0, more preferably at least about 1 .25, more preferably at 
least about 1.5 and most preferably at least about 2.0 indicates the action of positive 
selection and is considered to be an evolutionarily significant change. 

Strictly speaking, only K A /K S ratios greater than 1 .0 are indicative of positive 
selection. It is commonly accepted that the ESTs in GenBank® and other public 
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databases often suffer from some degree of sequencing error, and even a few incorrect 
nucleotides can influence K A /K S scores. Thus, all pairwise comparisons that involve 
public ESTs must be undertaken with care. Due to the errors inherent in the publicly 
available databases, it is possible that these errors could depress a K A /K S ratio below 1.0. 
5 For this reason, K A /Ks ratios between 0.75 and 1 .0 should be examined carefully in order 
to determine whether or not a sequencing error has obscured evidence of positive 
selection. Such errors may be discovered through sequencing methods that are designed 
to be highly accurate. 

The term "positive evolutionarily significant change" means an evolutionarily 
1 0 significant change in a particular species that results in an adaptive change that is positive 
as compared to other related species. Examples of positive evolutionarily significant 
changes are changes that have resulted in enhanced cognitive abilities or enhanced or 
unique physiological conditions in humans and adaptive changes in chimpanzees that 
have resulted in the ability of the chimpanzees infected with HIV to be resistant to 
15 progression to full-blown AIDS. 

The term "enhanced breast development" refers to the enlarged breasts observed 
in humans relative to non-human primates. The enlarged human breast has increased 
adipose, duct and/or gland tissue relative to other primates, and develops prior to first 
pregnancy and lactation. 
20 The term "resistant" means that an organism, such as a chimpanzee, exhibits an 

ability to avoid, or diminish the extent of, a disease condition and/or development of the 
disease, preferably when compared to non-resistant organisms, typically humans. For 
example, a chimpanzee is resistant to certain impacts of HIV and other viral infections, 
and/or it does not develop the ultimate disease - AIDS. 
25 The term "susceptibility" means that an organism, such as a human, fails to avoid, 

or diminish the extent of, a disease condition and/or development of the disease 
condition, preferably when compared to an organism that is known to be resistant, such as 
a non-human primate, such as chimpanzee. For example, a human is susceptible to 
certain impacts of HIV and other viral infections and/or development of the ultimate 
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disease — AIDS. 

It is understood that resistance and susceptibility vary from individual to 
individual, and that, for purposes of this invention, these terms also apply to a group of 
individuals within a species, and comparisons of resistance and susceptibility generally 
5 refer to overall, average differences between species, although intra-specific comparisons 
may be used. 

The term "homologous" or "homologue" or "ortholog" is known and well 
understood in the art and refers to related sequences that share a common ancestor and is 
determined based on degree of sequence identity. These terms describe the relationship 

10 between a gene found in one species and the corresponding or equivalent gene in another 
species. For purposes of this invention homologous sequences are compared. 
"Homologous sequences" or "homologues" or "orthologs" are thought, believed, or 
known to be functionally related. A functional relationship may be indicated in any one 
of a number of ways, including, but not limited to, (a) degree of sequence identity; (b) 

15 same or similar biological function. Preferably, both (a) and (b) are indicated. The 

degree of sequence identity may vary, but is preferably at least 50% (when using standard 
sequence alignment programs known in the art), more preferably at least 60%, more 
preferably at least about 75%, more preferably at least about 85%. Homology can be 
determined using software programs readily available in the art, such as those discussed 

20 in Current Protocols in Molecular Biology (F.M. Ausubel et aL, eds., 1987) Supplement 
30, section 7.718, Table 7.71 . Preferred alignment programs are Mac Vector (Oxford 
Molecular Ltd, Oxford, U.K.) and ALIGN Plus (Scientific and Educational Software, 
Pennsylvania). Another preferred alignment program is Sequencher (Gene Codes, Ann 
Arbor, Michigan), using default parameters. 

25 The term "nucleotide change" refers to nucleotide substitution, deletion, and/or 

insertion, as is well understood in the art. 

The term "human protein-coding nucleotide sequence" which is "associated with 
susceptibility to AIDS" as used herein refers to a human nucleotide sequence that encodes 
a protein that is associated with HIV dissemination (within the organism, i.e., intra- 
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organism infectivity), propagation and/or development of AIDS. Due to the extensive 
research in the mechanisms underlying progression from HIV infection to the 
development of AIDS, a number of candidate human genes are believed or known to be 
associated with one or more of these phenomena. A polynucleotide (including any 

5 polypeptide encoded therein) sequence associated with susceptibility to AIDS is one 
which is either known or implicated to play a role in HIV dissemination, replication, 
and/or subsequent progression to full-blown AIDS. Examples of such candidate genes 
are provided below. 

"AIDS resistant" means that an organism, such as a chimpanzee, exhibits an 

10 ability to avoid, or diminish the extent of, the result of HIV infection (such as propagation 
and dissemination) and/or development of AIDS, preferably when compared to AIDS- 
susceptible humans. 

"Susceptibility" to AIDS means that an organism, such as a human, fails to avoid, 
or diminish the extent of, the result of HIV infection (such as propagation and 

1 5 dissemination) and/or development of AIDS, preferably when compared to an organism 
that is known to be AIDS resistant, such as a non-human primate, such as chimpanzee. 

The term "brain protein-coding nucleotide sequence" as used herein refers to a 
nucleotide sequence expressed in the brain that encodes a protein. One example of the 
"brain protein-coding nucleotide sequence" is a brain cDNA sequence. 

20 As used herein, the term "brain functions unique or enhanced in humans" or 

"unique functional capabilities of the human brain" or "brain functional capability that is 
unique or enhanced in humans" refers to any brain function, either in kind or in degree, 
that is identified and/or observed to be enhanced in humans compared to other non- 
human primates. Such brain functions include, but are not limited to high capacity 

25 information processing, storage and retrieval capabilities, creativity, memory, language 
abilities, brain-mediated emotional response, locomotion, pain/pleasure sensation, 
olfaction, and temperament. 

"Housekeeping genes" is a term well understood in the art and means those genes 
associated with general cell function, including but not limited to growth, division, stasis, 
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metabolism, and/or death. "Housekeeping" genes generally perform functions found in 
more than one cell type. In contrast, cell-specific genes generally perform functions in a 
particular cell type (such as neurons) and/or class (such as neural cells). 

The term "agent", as used herein, means a biological or chemical compound such 
5 as a simple or complex organic or inorganic molecule, a peptide, a protein or an 

oligonucleotide. A vast array of compounds can be synthesized, for example oligomers, 
such as oligopeptides and oligonucleotides, and synthetic organic and inorganic 
compounds based on various core structures, and these are also included in the term 
"agent". In addition, various natural sources can provide compounds for screening, such 
1 0 as plant or animal extracts, and the like. Compounds can be tested singly or in 
combination with one another. 

The term "to modulate function" of a polynucleotide or a polypeptide means that 
the function of the polynucleotide or polypeptide is altered when compared to not adding 
an agent. Modulation may occur on any level that affects function. A polynucleotide or 
1 5 polypeptide function may be direct or indirect, and measured directly or indirectly. 

A "function of a polynucleotide" includes, but is not limited to, replication; 
translation; and expression pattern(s). A polynucleotide function also includes functions 
associated with a polypeptide encoded within the polynucleotide. For example, an agent 
which acts on a polynucleotide and affects protein expression, conformation, folding (or 
20 other physical characteristics), binding to other moieties (such as ligands), activity (or 
other functional characteristics), regulation and/or other aspects of protein structure or 
function is considered to have modulated polynucleotide function. 

A "function of a polypeptide" includes, but is not limited to, conformation, folding 
(or other physical characteristics), binding to other moieties (such as ligands), activity (or 
25 other functional characteristics), and/or other aspects of protein structure or functions. 
For example, an agent that acts on a polypeptide and affects its conformation, folding (or 
other physical characteristics), binding to other moieties (such as ligands), activity (or 
other functional characteristics), and/or other aspects of protein structure or functions is 
considered to have modulated polypeptide function. The ways that an effective agent can 
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act to modulate the function of a polypeptide include, but are not limited to 1) changing 
the conformation, folding or other physical characteristics; 2) changing the binding 
strength to its natural ligand or changing the specificity of binding to ligands; and 3) 
altering the activity of the polypeptide. 
5 The terms "modulate susceptibility to development of AIDS" and "modulate 

resistance to development of AIDS", as used herein, include modulating intra-organism 
cell-to-cell transmission or infectivity of HIV. The terms further include reducing 
susceptibility to development of AIDS and/or cell-to-cell transmission or infectivity of 
HIV. The terms further include increasing resistance to development of AIDS and/or 
1 0 cell-to-cell transmission or infectivity of HIV. One means of assessing whether an agent 
is one that modulates susceptibility or resistance to development of AIDS is to determine 
whether at least one index of HIV susceptibility is affected, using a cell-based system as 
described herein, as compared with an appropriate control. Indicia of HIV susceptibility 
include, but are not limited to, cell-to-cell transmission of the virus, as measured by total 
1 5 number of cells infected with HIV and syncytia formation. 

The term "target site" means a location in a polypeptide which can be one or more 
amino acids and/or is a part of a structural and/or functional motif, e.g., a binding site, a 
dimerization domain, or a catalytic active site. It also includes a location in a 
polynucleotide where there is one or more non-synonymous nucleotide changes in a 
20 protein coding region, or may also refer to a regulatory region of a positively selected 

gene. Target sites may be a useful for direct or indirect interaction with an agent, such as 
a therapeutic agent. 

The term "molecular difference" includes any structural and/or functional 
difference. Methods to detect such differences, as well as examples of such differences, 
25 are described herein. 

A "functional effect" is a term well known in the art, and means any effect which 
is exhibited on any level of activity, whether direct or indirect. 

General Procedures Known in the Art 
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For the purposes of this invention, the source of the human and non-human 
polynucleotide can be any suitable source, e.g., genomic sequences or cDNA sequences. 
Preferably, cDNA sequences from human and a non-human primate are compared. 
Human protein-coding sequences can be obtained from public databases such as the 

5 Genome Sequence Data Bank and GenBank. These databases serve as repositories of the 
molecular sequence data generated by ongoing research efforts. Alternatively, human 
protein-coding sequences may be obtained from, for example, sequencing of cDNA 
reverse transcribed from mRNA expressed in human cells, or after PCR amplification, 
according to methods well known in the art. Alternatively, human genomic sequences 

10 may be used for sequence comparison. Human genomic sequences can be obtained from 
public databases or from a sequencing of commercially available human genomic DNA 
libraries or from genomic DNA, after PCR. 

The non-human primate protein-coding sequences can be obtained by, for 
example, sequencing cDNA clones that are randomly selected from a non-human primate 

1 5 cDN A library. The non-human primate cDNA library can be constructed from total 
mRNA expressed in a primate cell using standard techniques in the art. In some 
embodiments, the cDNA is prepared from mRNA obtained from a tissue at a determined 
developmental stage, or a tissue obtained after the primate has been subjected to certain 
environmental conditions. cDNA libraries used for the sequence comparison of the 

20 present invention can be constructed using conventional cDNA library construction 

techniques that are explained fully in the literature of the art. Total mRNAs are used as 
templates to reverse-transcribe cDNAs. Transcribed cDNAs are subcloned into 
appropriate vectors to establish a cDNA library. The established cDNA library can be 
maximized for full-length cDNA contents, although less than full-length cDNAs maybe 

25 used. Furthermore, the sequence frequency can be normalized according to, for example, 
Bonaldo et al (1996) Genome Research 6:791-806. cDNA clones randomly selected 
from the constructed cDNA library can be sequenced using standard automated 
sequencing techniques. Preferably, full-length cDNA clones are used for sequencing. 
Either the entire or a large portion of cDNA clones from a cDNA library may be 
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sequenced, although it is also possible to practice some embodiments of the invention by 
sequencing as little as a single cDNA, or several cDNA clones. 

In one preferred embodiment of the present invention, non-human primate cDNA 
clones to be sequenced can be pre-selected according to their expression specificity. In 
5 order to select cDNAs corresponding to active genes that are specifically expressed, the 
cDNAs can be subject to subtraction hybridization using mRNAs obtained from other 
organs, tissues or cells of the same animal. Under certain hybridization conditions with 
appropriate stringency and concentration, those cDNAs that hybridize with non-tissue 
specific mRNAs and thus likely represent "housekeeping" genes will be excluded from 
1 0 the cDNA pool. Accordingly, remaining cDNAs to be sequenced are more likely to be 
associated with tissue-specific functions. For the purpose of subtraction hybridization, 
non-tissue-specific mRNAs can be obtained from one organ, or preferably from a 
combination of different organs and cells. The amount of non-tissue-specific mRNAs are 
maximized to saturate the tissue-specific cDNAs. 
1 5 Alternatively, information from online public databases can be used to select or 

give priority to cDNAs that are more likely to be associated with specific functions. For 
example, the non-human primate cDNA candidates for sequencing can be selected by 
PCR using primers designed from candidate human cDNA sequence. Candidate human 
cDNA sequences are, for example, those that are only found in a specific tissue, such as 
20 brain or breast, or that correspond to genes likely to be important in the specific function, 
such as brain function or breast tissue adipose or glandular development. Such human 
tissue-specific cDNA sequences can be obtained by searching online human sequence 
databases such as GenBank, in which information with respect to the expression profile 
and/or biological activity for cDNA sequences are specified. 
25 Sequences of non-human primate (for example, from an AIDS-resistant non- 

human primate) homologue(s) to a known human gene may be obtained using methods 
standard in the art, such as from public databases such as GenBank or PCR methods 
(using, for example, GeneAmp PCR System 9700 thermocyclers (Applied Biosystems, 
Inc.)). For example non-human primate cDNA candidates for sequencing can be selected 
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by PCR using primers designed from candidate human cDNA sequences. For PCR, 
primers may be made from the human sequences using standard methods in the art, 
including publicly available primer design programs such as PRIMER® (Whitehead 
Institute). The sequence amplified may then be sequenced using standard methods and 
5 equipment in the art, such as automated sequencers (Applied Biosystems, Inc.). 

GENERAL METHODS OF THE INVENTION 

The general method of the invention is as follows. Briefly, nucleotide sequences 
are obtained from a human source and a non-human source. The human and non-human 

10 nucleotide sequences are compared to one another to identify sequences that are 

homologous. The homologous sequences are analyzed to identify those that have nucleic 
acid sequence differences between the two species. Then molecular evolution analysis is 
conducted to evaluate quantitatively and qualitatively the evolutionary significance of the 
differences. For genes that have been positively selected between two species, e.g., 

1 5 human and chimp, it is useful to determine whether the difference occurs in other non- 
human primates. Next, the sequence is characterized in terms of molecular/genetic 
identity and biological function. Finally, the information can be used to identify agents 
useful in diagnosis and treatment of human medically or commercially relevant 
conditions. 

20 The general methods of the invention entail comparing human protein-coding 

nucleotide sequences to protein-coding nucleotide sequences of a non-human, preferably 
a primate, and most preferably a chimpanzee. Examples of other non-human primates are 
bonobo, gorilla, orangutan, gibbon, Old World monkeys, and New World monkeys. A 
phylogenetic tree for primates within the hominoid group is depicted in FIG. 1. 

25 Bioinformatics is applied to the comparison and sequences are selected that contain a 
nucleotide change or changes that is/are evolutionarily significant change(s). The 
invention enables the identification of genes that have evolved to confer some 
evolutionary advantage and the identification of the specific evolved changes. 
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Protein-coding sequences of human and another non-human primate are compared 
to identify homologous sequences. Protein-coding sequences known to or suspected of 
having a specific biological function may serve as the starting point for the comparison. 
Any appropriate mechanism for completing this comparison is contemplated by this 

5 invention. Alignment may be performed manually or by software (examples of suitable 
alignment programs are known in the art). Preferably, protein-coding sequences from a 
non-human primate are compared to human sequences via database searches, e.g., 
BLAST searches. The high scoring "hits," i.e., sequences that show a significant 
similarity after BLAST analysis, will be retrieved and analyzed. Sequences showing a 

1 0 significant similarity can be those having at least about 60%, at least about 75%, at least 
about 80%, at least about 85%, or at least about 90% sequence identity. Preferably, 
sequences showing greater than about 80% identity are further analyzed. The 
homologous sequences identified via database searching can be aligned in their entirety 
using sequence alignment methods and programs that are known and available in the art, 

1 5 such as the commonly used simple alignment program CLUSTAL V by Higgins et al 
(1992) GIB/OS 8:189-191. 

Alternatively, the sequencing and homologous comparison of protein-coding 
sequences between human and a non-human primate may be performed simultaneously 
by using the newly developed sequencing chip technology. See, for example, Rava et al 

20 US Patent 5,545,531. 

The aligned protein-coding sequences of human and another non-human primate 
are analyzed to identify nucleotide sequence differences at particular sites. Again, any 
suitable method for achieving this analysis is contemplated by this invention. If there are 
no nucleotide sequence differences, the non-human primate protein coding sequence is 

25 not usually further analyzed. The detected sequence changes are generally, and 

preferably, initially checked for accuracy. Preferably, the initial checking comprises 
performing one or more of the following steps, any and all of which are known in the art: 
(a) finding the points where there are changes between the non-human primate and 
human sequences; (b) checking the sequence fluorogram (chromatogram) to determine if 



the bases that appear unique to non-human primate correspond to strong, clear signals 
specific for the called base; (c) checking the human hits to see if there is more than one 
human sequence that corresponds to a sequence change. Multiple human sequence 
entries for the same gene that have the same nucleotide at a position where there is a 
different nucleotide in a non-human primate sequence provides independent support that 
the human sequence is accurate, and that the change is significant. Such changes are 
examined using public database information and the genetic code to determine whether 
these nucleotide sequence changes result in a change in the amino acid sequence of the 
encoded protein. As the definition of "nucleotide change" makes clear, the present 
invention encompasses at least one nucleotide change, either a substitution, a deletion or 
an insertion, in a human protein-coding polynucleotide sequence as compared to 
corresponding sequence from a non-human primate. Preferably, the change is a 
nucleotide substitution. More preferably, more than one substitution is present in the 
identified human sequence and is subjected to molecular evolution analysis. 

Any of several different molecular evolution analyses or K A /K s -type methods can 
be employed to evaluate quantitatively and qualitatively the evolutionary significance of 
the identified nucleotide changes between human gene sequences and that of a non- 
human primate. Kreitman and Akashi (1995) Annu. Rev. Ecol Syst. 26:403-422; Li, 
Molecular Evolution, Sinauer Associates, Sunderland, MA, 1997. For example, positive 
selection on proteins (i.e., molecular-level adaptive evolution) can be detected in protein- 
coding genes by pairwise comparisons of the ratios of nonsynonymous nucleotide 
substitutions per nonsynonymous site (K A ) to synonymous substitutions per synonymous 
site (K s ) (Li et al, 1985; Li, 1993). Any comparison of K A and K s may be used, although 
it is particularly convenient and most effective to compare these two variables as a ratio. 
Sequences are identified by exhibiting a statistically significant difference between K A 
and Ks using standard statistical methods. 

Preferably, the K A /K S analysis by Li et al is used to carry out the present 
invention, although other analysis programs that can detect positively selected genes 
between species can also be used. Li et al (1985) Mol Biol Evol 2:150-174; Li (1993); 
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see also J. Mol. Evol. 36:96-99; Messier and Stewart (1997) Nature 385:151-154; Nei 
(1987) Molecular Evolutionary Genetics (New York, Columbia University Press). The 
K A /K S method, which comprises a comparison of the rate of non-synonymous 
substitutions per non-synonymous site with the rate of synonymous substitutions per 
5 synonymous site between homologous protein-coding region of genes in terms of a ratio, 
is used to identify sequence substitutions that may be driven by adaptive selections as 
opposed to neutral selections during evolution. A synonymous ("silent") substitution is 
one that, owing to the degeneracy of the genetic code, makes no change to the amino acid 
sequence encoded; a non-synonymous substitution results in an amino acid replacement. 
1 0 The extent of each type of change can be estimated as K A and K s , respectively, the 
numbers of synonymous substitutions per synonymous site and non-synonymous 
substitutions per non-synonymous site. Calculations of K A /Ks may be performed 
manually or by using software. An example of a suitable program is MEGA (Molecular 
Genetics Institute, Pennsylvania State University). 
1 5 For the purpose of estimating K A and K s , either complete or partial human 

protein-coding sequences are used to calculate total numbers of synonymous and non- 
synonymous substitutions, as well as non-synonymous and synonymous sites. The length 
of the polynucleotide sequence analyzed can be any appropriate length. Preferably, the 
entire coding sequence is compared, in order to determine any and all significant changes. 
20 Publicly available computer programs, such as Li93 (Li (1993) J. Mol. Evol. 36:96-99) or 
INA, can be used to calculate the K A and K s values for all pairwise comparisons. This 
analysis can be further adapted to examine sequences in a "sliding window" fashion such 
that small numbers of important changes are not masked by the whole sequence. "Sliding 
window" refers to examination of consecutive, overlapping subsections of the gene (the 
25 subsections can be of any length). 

The comparison of non-synonymous and synonymous substitution rates is 
represented by the K A /K s ratio. K A /K s has been shown to be a reflection of the degree to 
which adaptive evolution has been at work in the sequence under study. Full length or 
partial segments of a coding sequence can be used for the K A /K s analysis. The higher the 
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K A /K s ratio, the more likely that a sequence has undergone adaptive evolution and the 
non-synonymous substitutions are evolutionarily significant. See, for example, Messier 
and Stewart (1997). Preferably, the K A /K S ratio is at least about 0.75, more preferably at 
least about 1.0, more preferably at least about 1.25, more preferably at least about 1.50, or 
5 more preferably at least about 2.00. Preferably, statistical analysis is performed on all 
elevated K A /K S ratios, including, but not limited to, standard methods such as Student's t- 
test and likelihood ratio tests described by Yang (1998) Mol. Biol Evol. 37:441-456. 

K A /K S ratios significantly greater than unity strongly suggest that positive 
selection has fixed greater numbers of amino acid replacements than can be expected as a 
1 0 result of chance alone, and is in contrast to the commonly observed pattern in which the 
ratio is less than or equal to one. Nei (1987); Hughes and Hei (1988) Nature 335:167- 
170; Messier and Stewart (1994) Current Biol. 4:911-913; Kreitman and Akashi (1995) 
Ann. Rev. Ecol. Syst. 26:403-422; Messier and Stewart (1997). Ratios less than one 
generally signify the role of negative, or purifying selection: there is strong pressure on 
1 5 the primary structure of functional, effective proteins to remain unchanged. 

All methods for calculating K A /K S ratios are based on a pairwise comparison of 
the number of nonsynonymous substitutions per nonsynonymous site to the number of 
synonymous substitutions per synonymous site for the protein-coding regions of 
homologous genes from related species. Each method implements different corrections 
20 for estimating "multiple hits" (i.e., more than one nucleotide substitution at the same site). 
Each method also uses different models for how DNA sequences change over 
evolutionary time. Thus, preferably, a combination of results from different algorithms is 
used to increase the level of sensitivity for detection of positively-selected genes and 
confidence in the result. 
25 Preferably, K A /K S ratios should be calculated for orthologous gene pairs, as 

opposed to paralogous gene pairs (i.e., a gene which results from speciation, as opposed 
to a gene that is the result of gene duplication) Messier and Stewart (1997). This 
distinction may be made by performing additional comparisons with other non-human 
primates, such as gorilla and orangutan, which allows for phylogenetic tree-building. 
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Orthologous genes when used in tree-building will yield the known "species tree", i.e., 
will produce a tree that recovers the known biological tree. In contrast, paralogous genes 
will yield trees which will violate the known biological tree. 

It is understood that the methods described herein could lead to the identification 
5 of human polynucleotide sequences that are functionally related to human protein-coding 
sequences. Such sequences may include, but are not limited to, non-coding sequences or 
coding sequences that do not encode human proteins. These related sequences can be, for 
example, physically adjacent to the human protein-coding sequences in the human 
genome, such as introns or 5'- and 3'- flanking sequences (including control elements 
10 such as promoters and enhancers). These related sequences may be obtained via 
searching a public human genome database such as GenBank or, alternatively, by 
screening and sequencing a human genomic library with a protein-coding sequence as 
probe. Methods and techniques for obtaining non-coding sequences using related coding 
sequence are well known to one skilled in the art. 
1 5 The evolutionarily significant nucleotide changes, which are detected by 

molecular evolution analysis such as the K A /K S analysis, can be further assessed for their 
unique occurrence in humans (or the non-human primate) or the extent to which these 
changes are unique in humans (or the non-human primate). For example, the identified 
changes can be tested for presence/absence in other non-human primate sequences. The 
20 sequences with at least one evolutionarily significant change between human and one 

non-human primate can be used as primers for PCR analysis of other non-human primate 
protein-coding sequences, and resulting polynucleotides are sequenced to see whether the 
same change is present in other non-human primates. These comparisons allow further 
discrimination as to whether the adaptive evolutionary changes are unique to the human 
25 lineage as compared to other non-human primates or whether the adaptive change is 
unique to the non-human primates (i.e., chimpanzee) as compared to humans and other 
non-human primates. A nucleotide change that is detected in human but not other 
primates more likely represents a human adaptive evolutionary change. Alternatively, a 
nucleotide change that is detected in a non-human primate (i.e., chimpanzee) that is not 
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detected in humans or other non-human primates likely represents a chimpanzee adaptive 
evolutionary change. Other non-human primates used for comparison can be selected 
based on their phylogenetic relationships with human. Closely related primates can be 
those within the hominoid sublineage, such as chimpanzee, bonobo, gorilla, and 
5 orangutan. Non-human primates can also be those that are outside the hominoid group 
and thus not so closely related to human, such as the Old World monkeys and New World 
monkeys. Statistical significance of such comparisons may be determined using 
established available programs, e.g., Mest as used by Messier and Stewart (1997) Nature 
385: 1 5 1-1 54. Those genes showing statistically high K A /K S ratios are very likely to have 
10 undergone adaptive evolution. 

Sequences with significant changes can be used as probes in genomes from 
different human populations to see whether the sequence changes are shared by more than 
one human population. Gene sequences from different human populations can be 
obtained from databases made available by, for example, the Human Genome Project, the 
15 human genome diversity project or, alternatively, from direct sequencing of PCR- 

amplified DNA from a number of unrelated, diverse human populations. The presence of 
the identified changes in different human populations would further indicate the 
evolutionary significance of the changes. Chimpanzee sequences with significant 
changes can be obtained and evaluated using similar methods to determine whether the 
20 sequence changes are shared among many chimpanzees. 

Sequences with significant changes between species can be further characterized 
in terms of their molecular/genetic identities and biological functions, using methods and 
techniques known to those of ordinary skill in the art. For example, the sequences can be 
located genetically and physically within the human genome using publicly available bio- 
25 informatics programs. The newly identified significant changes within the nucleotide 
sequence may suggest a potential role of the gene in human evolution and a potential 
association with human-unique functional capabilities. The putative gene with the 
identified sequences may be further characterized by, for example, homologue searching. 
Shared homology of the putative gene with a known gene may indicate a similar 
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biological role or function. Another exemplary method of characterizing a putative gene 
sequence is on the basis of known sequence motifs. Certain sequence patterns are known 
to code for regions of proteins having specific biological characteristics such as signal 
sequences, DNA binding domains, or transmembrane domains. 

The identified human sequences with significant changes can also be further 
evaluated by looking at where the gene is expressed in terms of tissue- or cell type- 
specificity. For example, the identified coding sequences can be used as probes to 
perform in situ mRNA hybridization that will reveal the expression patterns of the 
sequences. Genes that are expressed in certain tissues may be better candidates as being 
associated with important human functions associated with that tissue, for example brain 
tissue. The timing of the gene expression during each stage of human development can 
also be determined. 

As another exemplary method of sequence characterization, the functional roles of 
the identified nucleotide sequences with significant changes can be assessed by 
conducting functional assays for different alleles of an identified gene in a model system, 
such as yeast, nematode, Drosophila, and mouse. Model systems may be cell-based or in 
vivo, such as transgenic animals. Preferably, the transgenic mouse system is used. 
Methods of making cell-based systems and/or transgenic animal systems are known in the 
art and need not be described in detail herein. 

As another exemplary method of sequence characterization, the use of computer 
programs allows modeling and visualizing the three-dimensional structure of the 
homologous proteins from human and chimpanzee. Specific, exact knowledge of which 
amino acids have been replaced in the chimpanzee protein(s) allows detection of 
structural changes that may be associated with functional differences. Thus, use of 
modeling techniques is closely associated with identification of functional roles discussed 
in the previous paragraph. The use of individual or combinations of these techniques 
constitutes part of the present invention. For example, chimpanzee ICAM-3 contains a 
glutamine residue (Q101) at the site in which human ICAM-3 contains a proline (P101). 
The human protein is known to bend sharply at this point. Replacement of the proline by 
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glutamine in the chimpanzee protein is likely to result in a much less sharp bend at this 
point. This has clear implications for packaging of the IC AM-3 chimpanzee protein into 
HIV virions. 

The sequences identified by the methods described herein have significant uses in 
5 diagnosis and treatment of medically or commercially relevant human conditions. 

Accordingly, the present invention provides methods for identifying agents that are useful 
in modulating human-unique or human-enhanced functional capabilities and/or correcting 
defects in these capabilities using these sequences. These methods employ, for example, 
screening techniques known in the art, such as in vitro systems, cell-based expression 
10 systems and transgenic animal systems. The approach provided by the present invention 
not only identifies rapidly evolved genes, but indicates modulations that can be made to 
the protein that may not be too toxic because they exist in another species. 
Screening Methods 

The present invention also provides screening methods using the polynucleotides 
1 5 and polypeptides identified and characterized using the above-described methods. These 
screening methods are useful for identifying agents which may modulate the function(s) 
of the polynucleotides or polypeptides in a manner that would be useful for a human 
treatment. Generally, the methods entail contacting at least one agent to be tested with 
either a cell that has been transfected with a polynucleotide sequence identified by the 
20 methods described above, or a preparation of the polypeptide encoded by such 

polynucleotide sequence, wherein an agent is identified by its ability to modulate function 
of either the polynucleotide sequence or the polypeptide. 

As used herein, the term "agent" means a biological or chemical compound such 
as a simple or complex organic or inorganic molecule, a peptide, a protein or an 
25 oligonucleotide. A vast array of compounds can be synthesized, for example oligomers, 
such as oligopeptides and oligonucleotides, and synthetic organic and inorganic 
compounds based on various core structures, and these are also included in the term 
"agent". In addition, various natural sources can provide compounds for screening, such 
as plant or animal extracts, and the like. Compounds can be tested singly or in 
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combination with one another. 

To "modulate function" of a polynucleotide or a polypeptide means that the 
function of the polynucleotide or polypeptide is altered when compared to not adding an 
agent. Modulation may occur on any level that affects function. A polynucleotide or 
5 polypeptide function may be direct or indirect, and measured directly or indirectly. A 
"function" of a polynucleotide includes, but is not limited to, replication, translation, and 
expression pattern(s). A polynucleotide function also includes functions associated with 
a polypeptide encoded within the polynucleotide. For example, an agent which acts on a 
polynucleotide and affects protein expression, conformation, folding (or other physical 
10 characteristics), binding to other moieties (such as ligands), activity (or other functional 
characteristics), regulation and/or other aspects of protein structure or function is 
considered to have modulated polynucleotide function. The ways that an effective agent 
can act to modulate the expression of a polynucleotide include, but are not limited to 1) 
modifying binding of a transcription factor to a transcription factor responsive element in 
15 the polynucleotide; 2) modifying the interaction between two transcription factors 

necessary for expression of the polynucleotide; 3) altering the ability of a transcription 
factor necessary for expression of the polynucleotide to enter the nucleus; 4) inhibiting 
the activation of a transcription factor involved in transcription of the polynucleotide; 5) 
modifying a cell-surface receptor which normally interacts with a ligand and whose 
20 binding of the ligand results in expression of the polynucleotide; 6) inhibiting the 

inactivation of a component of the signal transduction cascade that leads to expression of 
the polynucleotide; and 7) enhancing the activation of a transcription factor involved in 
transcription of the polynucleotide. 

A "function" of a polypeptide includes, but is not limited to, conformation, folding 
25 (or other physical characteristics), binding to other moieties (such as ligands), activity (or 
other functional characteristics), and/or other aspects of protein structure or functions. 
For example, an agent that acts on a polypeptide and affects its conformation, folding (or 
other physical characteristics), binding to other moieties (such as ligands), activity (or 
other functional characteristics), and/or other aspects of protein structure or functions is 
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considered to have modulated polypeptide function. The ways that an effective agent can 
act to modulate the function of a polypeptide include, but are not limited to 1) changing 
the conformation, folding or other physical characteristics; 2) changing the binding 
strength to its natural ligand or changing the specificity of binding to ligands; and 3) 
altering the activity of the polypeptide. 

A "function" of a polynucleotide includes its expression, i.e., transcription and/or 
translation. It can also include (without limitation) its conformation, folding and binding 
to other moieties. 

Generally, the choice of agents to be screened is governed by several parameters, 
such as the particular polynucleotide or polypeptide target, its perceived function, its 
three-dimensional structure (if known or surmised), and other aspects of rational drug 
design. Techniques of combinatorial chemistry can also be used to generate numerous 
permutations of candidates. Those of skill in the art can devise and/or obtain suitable 
agents for testing. 

The in vivo screening assays described herein may have several advantages over 
conventional drug screening assays: 1) if an agent must enter a cell to achieve a desired 
therapeutic effect, an in vivo assay can give an indication as to whether the agent can 
enter a cell; 2) an in vivo screening assay can identify agents that, in the state in which 
they are added to the assay system are ineffective to elicit at least one characteristic which 
is associated with modulation of polynucleotide or polypeptide function, but that are 
modified by cellular components once inside a cell in such a way that they become 
effective agents; 3) most importantly, an in vivo assay system allows identification of 
agents affecting any component of a pathway that ultimately results in characteristics that 
are associated with polynucleotide or polypeptide function. 

In general, screening can be performed by adding an agent to a sample of 
appropriate cells which have been transfected with a polynucleotide identified using the 
methods of the present invention, and monitoring the effect, i.e., modulation of a function 
of the polynucleotide or the polypeptide encoded within the polynucleotide. The 
experiment preferably includes a control sample which does not receive the candidate 
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agent. The treated and untreated cells are then compared by any suitable phenotypic 
criteria, including but not limited to microscopic analysis, viability testing, ability to 
replicate, histological examination, the level of a particular RNA or polypeptide 
associated with the cells, the level of enzymatic activity expressed by the cells or cell 
5 lysates, the interactions of the cells when exposed to infectious agents, such as HIV, and 
the ability of the cells to interact with other cells or compounds. For example, the 
transfected cells can be exposed to the agent to be tested and, before, during, or after 
treatment with the agent, the cells can be infected with a virus, such as HIV, and tested 
for any indication of susceptibility of the cells to viral infection, including, for example, 
1 0 susceptibility of the cells to cell-to-cell viral infection, replication of the virus, production 
of a viral protein, and/or syncytia formation following infection with the virus. 
Differences between treated and untreated cells indicate effects attributable to the 
candidate agent. Optimally, the agent has a greater effect on experimental cells than on 
control cells. Appropriate host cells include, but are not limited to, eukaryotic cells, 
1 5 preferably mammalian cells. The choice of cell will at least partially depend on the 
nature of the assay contemplated. 

To test for agents that upregulate the expression of a polynucleotide, a suitable 
host cell transfected with a polynucleotide of interest, such that the polynucleotide is 
expressed (as used herein, expression includes transcription and/or translation) is 
20 contacted with an agent to be tested. An agent would be tested for its ability to result in 
increased expression of mRNA and/or polypeptide. Methods of making vectors and 
transfection are well known in the art. "Transfection" encompasses any method of 
introducing the exogenous sequence, including, for example, lipofection, transduction, 
infection or electroporation. The exogenous polynucleotide may be maintained as a non- 
25 integrated vector (such as a plasmid) or may be integrated into the host genome. 

To identify agents that specifically activate transcription, transcription regulatory 
regions could be linked to a reporter gene and the construct added to an appropriate host 
cell. As used herein, the term "reporter gene" means a gene that encodes a gene product 
that can be identified (i.e., a reporter protein). Reporter genes include, but are not limited 
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to, alkaline phosphatase, chloramphenicol acetyltransferase, P-galactosidase, luciferase 
and green fluorescence protein (GFP). Identification methods for the products of reporter 
genes include, but are not limited to, enzymatic assays and fluorimetric assays. Reporter 
genes and assays to detect their products are well known in the art and are described, for 

5 example in Ausubel et al. (1987) and periodic updates. Reporter genes, reporter gene 
assays, and reagent kits are also readily available from commercial sources. Examples of 
appropriate cells include, but are not limited to, fungal, yeast, mammalian, and other 
eukaryotic cells. A practitioner of ordinary skill will be well acquainted with techniques 
for transfecting eukaryotic cells, including the preparation of a suitable vector, such as a 

10 viral vector; conveying the vector into the cell, such as by electroporation; and selecting 
cells that have been transformed, such as by using a reporter or drug sensitivity element. 
The effect of an agent on transcription from the regulatory region in these constructs 
would be assessed through the activity of the reporter gene product. 

Besides the increase in expression under conditions in which it is normally 

1 5 repressed mentioned above, expression could be decreased when it would normally be 
maintained or increased. An agent could accomplish this through a decrease in 
transcription rate and the reporter gene system described above would be a means to 
assay for this. The host cells to assess such agents would need to be permissive for 
expression. 

20 Cells transcribing mRN A (from the polynucleotide of interest) could be used to 

identify agents that specifically modulate the half-life of mRN A and/or the translation of 
mRNA. Such cells would also be used to assess the effect of an agent on the processing 
and/or post-translational modification of the polypeptide. An agent could modulate the 
amount of polypeptide in a cell by modifying the turn-over (i.e., increase or decrease the 

25 half-life) of the polypeptide. The specificity of the agent with regard to the mRNA and 
polypeptide would be determined by examining the products in the absence of the agent 
and by examining the products of unrelated mRNAs and polypeptides. Methods to 
examine mRNA half-life, protein processing, and protein turn-over are well know to 
those skilled in the art. 
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In vivo screening methods could also be useful in the identification of agents that 
modulate polypeptide function through the interaction with the polypeptide directly. 
Such agents could block normal polypeptide-ligand interactions, if any, or could enhance 
or stabilize such interactions. Such agents could also alter a conformation of the 
5 polypeptide. The effect of the agent could be determined using immunoprecipitation 
reactions. Appropriate antibodies would be used to precipitate the polypeptide and any 
protein tightly associated with it. By comparing the polypeptides immunoprecipitated 
from treated cells and from untreated cells, an agent could be identified that would 
augment or inhibit polypeptide-ligand interactions, if any. Polypeptide-ligand 
10 interactions could also be assessed using cross-linking reagents that convert a close, but 
noncovalent interaction between polypeptides into a covalent interaction. Techniques to 
examine protein-protein interactions are well known to those skilled in the art. 
Techniques to assess protein conformation are also well known to those skilled in the art. 
It is also understood that screening methods can involve in vitro methods, such as 
1 5 cell-free transcription or translation systems. In those systems, transcription or translation 
is allowed to occur, and an agent is tested for its ability to modulate function. For an 
assay that determines whether an agent modulates the translation of mRNA or a 
polynucleotide, an in vitro transcription/translation system maybe used. These systems 
are available commercially and provide an in vitro means to produce mRNA 
20 corresponding to a polynucleotide sequence of interest. After mRNA is made, it can be 
translated in vitro and the translation products compared. Comparison of translation 
products between an in vitro expression system that does not contain any agent (negative 
control) with an in vitro expression system that does contain an agent indicates whether 
the agent is affecting translation. Comparison of translation products between control 
25 and test polynucleotides indicates whether the agent, if acting on this level, is selectively 
affecting translation (as opposed to affecting translation in a general, non-selective or 
non-specific fashion). The modulation of polypeptide function can be accomplished in 
many ways including, but not limited to, the in vivo and in vitro assays listed above as 
well as in in vitro assays using protein preparations. Polypeptides can be extracted and/or 
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purified from natural or recombinant sources to create protein preparations. An agent can 
be added to a sample of a protein preparation and the effect monitored; that is whether 
and how the agent acts on a polypeptide and affects its conformation, folding (or other 
physical characteristics), binding to other moieties (such as ligands), activity (or other 

5 functional characteristics), and/or other aspects of protein structure or functions is 
considered to have modulated polypeptide function. 

In an example for an assay for an agent that binds to a polypeptide encoded by a 
polynucleotide identified by the methods described herein, a polypeptide is first 
recombinantly expressed in a prokaryotic or eukaryotic expression system as a native or 

10 as a fusion protein in which a polypeptide (encoded by a polynucleotide identified as 
described above) is conjugated with a well-characterized epitope or protein. 
Recombinant polypeptide is then purified by, for instance, immunoprecipitation using 
appropriate antibodies or anti-epitope antibodies or by binding to immobilized ligand of 
the conjugate. An affinity column made of polypeptide or fusion protein is then used to 

15 screen a mixture of compounds which have been appropriately labeled. Suitable labels 
include, but are not limited to fluorochromes, radioisotopes, enzymes and 
chemiluminescent compounds. The unbound and bound compounds can be separated by 
washes using various conditions (e.g. high salt, detergent ) that are routinely employed by 
those skilled in the art. Non-specific binding to the affinity column can be minimized by 

20 pre-clearing the compound mixture using an affinity column containing merely the 

conjugate or the epitope. Similar methods can be used for screening for an agent(s) that 
competes for binding to polypeptides. In addition to affinity chromatography, there are 
other techniques such as measuring the change of melting temperature or the fluorescence 
anisotropy of a protein which will change upon binding another molecule. For example, 

25 a BIAcore assay using a sensor chip (supplied by Pharmacia Biosensor, Stitt et al (1995) 
Cell 80: 661-670) that is covalently coupled to polypeptide may be performed to 
determine the binding activity of different agents. 

It is also understood that the in vitro screening methods of this invention include 
structural, or rational, drug design, in which the amino acid sequence, three-dimensional 
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atomic structure or other property (or properties) of a polypeptide provides a basis for 
designing an agent which is expected to bind to a polypeptide. Generally, the design 
and/or choice of agents in this context is governed by several parameters, such as side-by- 
side comparison of the structures of a human and homologous non-human primate 
5 polypeptides, the perceived function of the polypeptide target, its three-dimensional 
structure (if known or surmised), and other aspects of rational drug design. Techniques 
of combinatorial chemistry can also be used to generate numerous permutations of 
candidate agents. 

Also contemplated in screening methods of the invention are transgenic animal 

10 systems, which are known in the art. 

The screening methods described above represent primary screens, designed to 
detect any agent that may exhibit activity that modulates the function of a polynucleotide 
or polypeptide. The skilled artisan will recognize that secondary tests will likely be 
necessary in order to evaluate an agent further. For example, a secondary screen may 

15 comprise testing the agent(s) in an infectivity assay using mice and other animal models 
(such as rat), which are known in the art. In addition, a cytotoxicity assay would be 
performed as a further corroboration that an agent which tested positive in a primary 
screen would be suitable for use in living organisms. Any assay for cytotoxicity would be 
suitable for this purpose, including, for example the MTT assay (Promega). 

20 The invention also includes agents identified by the screening methods described 

herein. 

Methods Useful for Identifying Positively Selected Non-Human Traits 

In one aspect of the invention, a non-human primate polynucleotide or 
25 polypeptide has undergone natural selection that resulted in a positive evolutionarily 
significant change (i.e., the non-human primate polynucleotide or polypeptide has a 
positive attribute not present in humans). In this aspect of the invention, the positively 
selected polynucleotide or polypeptide maybe associated with susceptibility or resistance 
to certain diseases or with other commercially relevant traits. Examples of this 
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embodiment include, but are not limited to, polynucleotides and polypeptides that have 
been positively selected in non-human primates, preferably chimpanzees, that may be 
associated with susceptibility or resistance to infectious diseases, cancer, or acne or may 
be associated with aesthetic conditions of interest to humans, such as hair growth or 

5 muscle mass. An example of this embodiment includes polynucleotides and polypeptides 
associated with the susceptibility or resistance to HIV progression to AIDS. The present 
invention can thus be useful in gaining insight into the molecular mechanisms that 
underlie resistance to HIV infection progressing to development of AIDS, providing 
information that can also be useful in discovering and/or designing agents such as drugs 

10 that prevent and/or delay development of AIDS. For example, CD59, which has been 
identified as a leukocyte and erythrocyte protein whose function is to protect these cells 
from the complement arm of the body's MAC (membrane attack complex) defense 
system (Meri et al (1996) Biochem. J. 616:923-935), has been found to be positively 
selected in the chimpanzee (see Example 16). It is believed that the CD59 found in 

15 chimpanzees confers a resistance to the progression of AIDS that is not found in humans. 
Thus, the positively selected chimpanzee CD59 can serve in the development of agents 
or drugs that are useful in arresting the progression of AIDS in humans, as is described in 
the Examples. 

Commercially relevant examples include, but are not limited to, polynucleotides 
20 and polypeptides that are positively selected in non-human primates that may be 
associated with aesthetic traits, such as hair growth, acne, or muscle mass. 

Accordingly, in one aspect, the invention provides methods for identifying a 
polynucleotide sequence encoding a polypeptide, wherein said polypeptide may be 
associated with a medically or commercially relevant positive evolutionarily significant 
25 change. The method comprises the steps of: (a) comparing human protein-coding 

nucleotide sequences to protein-coding nucleotide sequences of a non-human primate; 
and (b) selecting a non-human primate polynucleotide sequence that contains at least one 
nucleotide change as compared to corresponding sequence of the human, wherein said 
change is evolutionarily significant. The sequences identified by this method may be 
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further characterized and/or analyzed for their possible association with biologically or 
medically relevant functions unique or enhanced in non-human primates. 

Methods Useful for Identifying Positively Selected Human Traits 
5 This invention specifically provides methods for identifying human 

polynucleotide and polypeptide sequences that may be associated with unique or 
enhanced functional capabilities or traits of the human, for example, brain function or 
longer life span. More particularly, these methods identify those genetic sequences that 
may be associated with capabilities that are unique or enhanced in humans, including, but 

10 not limited to, brain functions such as high capacity information processing, storage and 
retrieval capabilities, creativity, and language abilities. Moreover, these methods identify 
those sequences that may be associated to other brain functional features with respect to 
which the human brain performs at enhanced levels as compared to other non-human 
primates; these differences may include brain-mediated emotional response, locomotion, 

15 pain/pleasure sensation, olfaction, temperament and longer life span. 

In this method, the general methods of the invention are applied as described 
above. Generally, the methods described herein entail (a) comparing human protein- 
coding polynucleotide sequences to that of a non-human primate; and (b) selecting those 
human protein-coding polynucleotide sequences having evolutionarily significant changes 

20 that may be associated with unique or enhanced functional capabilities of the human as 
compared to that of the non-human primate. 

In this embodiment, the human sequence includes the evolutionarily significant 
change (i.e., the human sequence differs from more than one non-human primate species 
sequence in a manner that suggests that such a change is in response to a selective 

25 pressure). The identity and function of the protein encoded by the gene that contains the 
evolutionarily significant change is characterized and a determination is made whether or 
not the protein can be involved in a unique or enhanced human function. If the protein is 
involved in a unique or enhanced human function, the information is used in a manner to 
identify agents that can supplement or otherwise modulate the unique or enhanced human 
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function. 

As a non-limiting example of the invention, identifying the genetic (i.e., 
nucleotide sequence) differences underlying the functional uniqueness of human brain 
may provide a basis for designing agents that can modulate human brain functions and/or 

5 help correct functional defects. These sequences could also be used in developing 
diagnostic reagents and/or biomedical research tools. The invention also provides 
methods for a large-scale comparison of human brain protein-coding sequences with 
those from a non-human primate. 

The identified human sequence changes can be used in establishing a database of 

10 candidate human genes that may be involved in human brain function. Candidates are 
ranked as to the likelihood that the gene is responsible for the unique or enhanced 
functional capabilities found in the human brain compared to chimpanzee or other non- 
human primates. Moreover, the database not only provides an ordered collection of 
candidate genes, it also provides the precise molecular sequence differences that exist 

1 5 between human and chimpanzee (and other non-human primates), and thus defines the 
changes that underlie the functional differences. This information can be useful in the 
identification of potential sites on the protein that may serve as useful targets for 
pharmaceutical agents. 

Accordingly, the present invention also provides methods for correlating an 

20 evolutionarily significant nucleotide change to a brain functional capability that is unique 
or enhanced in humans, comprising (a) identifying a human nucleotide sequence 
according to the methods described above; and (b) analyzing the functional effect of the 
presence or absence of the identified sequence in a model system. 

Further studies can be carried out to confirm putative function. For example, the 

25 putative function can be assayed in appropriate in vitro assays using transiently or stably 
transfected mammalian cells in culture, or using mammalian cells transfected with an 
antisense clone to inhibit expression of the identified polynucleotide to assess the effect 
of the absence of expression of its encoded polypeptide. Studies such as one-hybrid and 
two-hybrid studies can be conducted to determine, for example, what other 
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macromolecules the polypeptide interacts with. Transgenic nematodes or Drosophila can 
be used for various functional assays, including behavioral studies. The appropriate 
studies depend on the nature of the identified polynucleotide and the polypeptide encoded 
within the polynucleotide, and would be obvious to those skilled in the art. 

5 The present invention also provides polynucleotides and polypeptides identified 

by the methods of the present invention. In one embodiment, the present invention 
provides an isolated AATYK nucleotide sequence selected from the group consisting of 
nucleotides 2180-2329 of SEQ ED NO:14, nucleotides 2978-3478 of SEQ ID NO: 14, and 
nucleotides 3380-3988 of SEQ ID NO: 14; and an isolated nucleotide sequence having at 

1 0 least 85% homology to a nucleotide sequence of any of the preceding sequences. 

In another embodiment, the invention provides an isolated AATYK polypeptide 
selected from the group consisting of a polypeptide encoded by a nucleotide sequence 
selected from the group consisting of SEQ ID NO: 17 and SEQ ID NO: 18; wherein said 
encoding is based on the open reading frame (ORF) of SEQ ID NO: 14, and a polypeptide 

15 encoded by a nucleotide sequence having at least 85% homology to a nucleotide sequence 
selected from the group consisting of SEQ ID NO:17 and SEQ ID NO:18; wherein said 
encoding is based on the open reading frame of SEQ ID NO: 14. 

In a further embodiment, the present invention provides an isolated AATYK 
polypeptide selected from the group consisting of a polypeptide encoded by a nucleotide 

20 sequence selected from the group consisting of nucleotides 1-501 of SEQ ID NO: 17, 

nucleotides 1-150 of SEQ ID NO:17, nucleotides 100-249 of SEQ ID NO:17, nucleotides 
202-351 of SEQ ID NO:17, nucleotides 301-450 of SEQ ID NO:17, nucleotides 799-948 
of SEQ ID NO:17, nucleotides 901-1050 of SEQ ID NO:17, nucleotides 799-1299 of 
SEQ ID NO:17, and nucleotides 1201-1809 of SEQ ID NO:17; wherein said encoding is 

25 based on the open reading frame of SEQ ID NO: 14; and a polypeptide encoded by a 
nucleotide sequence having at least 85% homology to any of the preceding nucleotide 
sequences. 

In still another embodiment, the invention provides an isolated polypeptide 
selected from the group consisting of a polypeptide encoded by a nucleotide sequence 



selected from the group consisting of nucleotides 1-501 of SEQ ID NO:18, nucleotides 
799-1299 of SEQ ID NO: 18, and nucleotides 1201-1809 of SEQ ID NO:18; wherein said 
encoding is based on the open reading frame of SEQ ID NO: 14; and a polypeptide 
encoded by a nucleotide sequence having at least 85% homology to nucleotides 1-501 of 
5 SEQ ID NO:18, nucleotides 799-1299 of SEQ ID NO: 18, and nucleotides 1201-1809 of 
SEQK)NO:18. 

In another embodiment, the invention provides an isolated polynucleotide 
comprising SEQ ID NO: 17, wherein the coding capacity of the nucleic acid molecule is 
based on the open reading frame of SEQ ID NO: 14. In a preferred embodiment, the 

10 polynucleotide is a Pan troglodytes polynucleotide. 

In another embodiment, the invention provides an isolated polynucleotide 
comprising SEQ ID NO: 18, wherein the coding capacity of the nucleic acid molecule is 
based on the open reading frame of SEQ ID NO: 14. In a preferred embodiment, the 
polynucleotide is a Gorilla gorilla polynucleotide. 

15 In some embodiments, the polynucleotide or polypeptide having 85% homology 

to an isolated AATYK polynucleotide or polypeptide of the present invention is a 
homo log, which, when compared to a non-human primate, yields a Ka/K s ratio of at least 
0.75, at least 1.00, at least 1.25, at least 1.50, or at least 2.00. 

In other embodiments, the polynucleotide or polypeptide having 85% homology to 

20 an isolated AATYK polynucleotide or polypeptide of the present invention is a homolog 
which is capable of performing the function of the natural AATYK polynucleotide or 
polypeptide in a functional assay. Suitable assays for assessing the function of an 
ATTYK polynucleotide or polypeptide include a neuronal differentiation assay such as 
that described by Raghunath, et al., Brain Res Mol Brain Res. (2000) 77:151-62, or a 

25 tyrosine phosphorylation assay such as that described in Tomomura, et al, Oncogene 
(2001) 20(9):1022-32. The phrase "capable of performing the function of the natural 
AATYK polynucleotide or polypeptide in a functional assay" means that the 
polynucleotide or polypeptide has at least about 10% of the activity of the natural 
polynucleotide or polypeptide in the functional assay. In other preferred embodiments, 
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has at least about 20% of the activity of the natural polynucleotide or polypeptide in the 
functional assay. In other preferred embodiments, has at least about 30% of the activity 
of the natural polynucleotide or polypeptide in the functional assay. In other preferred 
embodiments, has at least about 40% of the activity of the natural polynucleotide or 

5 polypeptide in the functional assay. In other preferred embodiments, has at least about 
50% of the activity of the natural polynucleotide or polypeptide in the functional assay. 
In other preferred embodiments, the polynucleotide or polypeptide has at least about 60% 
of the activity of the natural polynucleotide or polypeptide in the functional assay. In 
more preferred embodiments, the polynucleotide or polypeptide has at least about 70% of 

1 0 the activity of the natural polynucleotide or polypeptide in the functional assay. In more 
preferred embodiments, the polynucleotide or polypeptide has at least about 80% of the 
activity of the natural polynucleotide or polypeptide in the functional assay. In more 
preferred embodiments, the polynucleotide or polypeptide has at least about 90% of the 
activity of the natural polynucleotide or polypeptide in the functional assay. 



Description of the AIDS Embodiment (an example of a positively selected non-human 
trait) 

The AIDS (Acquired Immune Deficiency Syndrome) epidemic has been estimated 
20 to threaten 30 million people world-wide (UNAIDS/WHO, 1998, "Report on the global 
HIV/ AIDS epidemic"). Well over a million people are infected in developed countries, 
and in parts of sub-Saharan Africa, 1 in 4 adults now carries the virus (UNAIDS/WHO, 
1998). Although efforts to develop vaccines are underway, near term prospects for 
successful vaccines are grim. Baiter and Cohen (1998) Science 281:159-160; Baltimore 
25 and Heilman (1998) Scientific Am. 279:98-103. Further complicating the development of 
therapeutics is the rapid mutation rate of HIV (the human immunodeficiency virus which 
is responsible for AIDS), which generates rapid changes in viral proteins. These changes 
ultimately allow the virus to escape current therapies, which target viral proteins. Dobkin 
(1998) Inf. Med. 15(3):159. Even drug cocktails which initially showed great promise are 
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subject to the emergence of drug-resistant mutants. Baiter and Cohen (1998); Dobkin 
(1998). Thus, there is still a serious need for development of therapies which delay or 
prevent progression of AIDS in HIV-infected individuals. Chun et al (1997) Proc. Natl 
Acad, ScL USA 94:13193-13197; Dobkin (1998). 
5 Human's closest relatives, chimpanzees (Pan troglodytes), have unexpectedly 

proven to be poor models for the study of the disease processes following infection with 
HIV-1. Novembre a/. (1997); J. Virol. 71(5):4086-4091. Once infected with HIV- 1, 
chimpanzees display resistance to progression of the disease. To date, only one 
chimpanzee individual is known to have developed full-blown AIDS, although more than 
10 100 captive chimpanzees have been infected. Novembre et al (1997); Villinger et al 
(1997) J. Med, Primatol 26(1-2): 1 1-18. Clearly, an understanding of the mechanism(s) 
that confer resistance to progression of the disease in chimpanzees may prove invaluable 
for efforts to develop therapeutic agents for HIV-infected humans. 

It is generally believed that wild chimpanzee populations harbored the HIV-1 
1 5 virus (perhaps for millennia) prior to its recent cross-species transmission to humans. 

Dube et al, (1994); Virology 202:379-389; Zhu and Ho (1995) Nature 374:503-504; Zhu 
et al (1998); Quinn (1994) Proc. Natl Acad. Sci USA 91:2407-2414. During this 
extended period, viral/host co-evolution has apparently resulted in accommodation, 
explaining chimpanzee resistance to ADDS progression. Burnet and White (1972); 
20 Natural History of Infectious Disease (Cambridge, Cambridge Univ. Press); Ewald 
(1991) Hum. Nat 2(i):l-30. All references cited herein are hereby incorporated by 
reference in their entirety. 

One aspect of this invention arises from the observations that (a) because 
chimpanzees (Pan troglodytes) have displayed resistance to development of AIDS 
25 although susceptible to HIV infection (Alter et al (1984) Science 226:549-552; Fultz et 
al (1986) 1 Virol. 58:1 16-124; Novembre et al (1997) 1 Virol 71(5):4086-4091), while 
humans are susceptible to developing this devastating disease, certain genes in 
chimpanzees may contribute to this resistance; and (b) it is possible to evaluate whether 
changes in human genes when compared to homologous genes from other species (such 
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as chimpanzee) are evolutionarily significant (z.e., indicating positive selective pressure). 
Thus, protein coding polynucleotides may contain sequence changes that are found in 
chimpanzees (as well as other AIDS-resistant primates) but not in humans, likely as a 
result of positive adaptive selection during evolution. Furthermore, such evolutionarily 

5 significant changes in polynucleotide and polypeptide sequences may be attributed to an 
AIDS-resistant non-human primate's (such as chimpanzee) ability to resist development 
of AIDS. The methods of this invention employ selective comparative analysis to 
identify candidate genes which may be associated with susceptibility or resistance to 
AIDS, which may provide new host targets for therapeutic intervention as well as specific 

10 information on the changes that evolved to confer resistance. Development of therapeutic 
approaches that involve host proteins (as opposed to viral proteins and/or mechanisms) 
may delay or even avoid the emergence of resistant viral mutants. The invention also 
provides screening methods using the sequences and structural differences identified. 
This invention provides methods for identifying human polynucleotide and 

1 5 polypeptide sequences that may be associated with susceptibility to post-infection 

development of AIDS. Conversely, the invention also provides methods for identifying 
polynucleotide and polypeptide sequences from an AIDS-resistant non-human primate 
(such as chimpanzee) that may be associated with resistance to development of AIDS. 
Identifying the genetic (i.e., nucleotide sequence) and the resulting protein structural and 

20 biochemical differences underlying susceptibility or resistance to development of ADDS 
will likely provide a basis for discovering and/or designing agents that can provide 
prevention and/or therapy for HIV infection progressing to AIDS. These differences 
could also be used in developing diagnostic reagents and/or biomedical research tools. 
For example, identification of proteins which confer resistance may allow development of 

25 diagnostic reagents or biomedical research tools based upon the disruption of the disease 
pathway of which the resistant protein plays a part. 

Generally, the methods described herein entail (a) comparing human protein- 
coding polynucleotide sequences to that of an AIDS resistant non-human primate (such as 
chimpanzee), wherein the human protein coding polynucleotide sequence is associated 
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with development of AIDS; and (b) selecting those human protein-coding polynucleotide 
sequences having evolutionarily significant changes that may be associated with 
susceptibility to development of AIDS. In another embodiment, the methods entail (a) 
comparing human protein-coding polynucleotide sequences to that of an AIDS-resistant 
5 non-human primate (such as chimpanzee), wherein the human protein coding 

polynucleotide sequence is associated with development of AIDS; and (b) selecting those 
non-human primate protein-coding polynucleotide sequences having evolutionarily 
significant changes that may be associated with resistance to development of AIDS. 
As is evident, the methods described herein can be applied to other infectious 

10 diseases. For example, the methods could be used in a situation in which a non-human 
primate is known or believed to have harbored the infectious disease for a significant 
period (i.e., a sufficient time to have allowed positive selection) and is resistant to 
development of the disease. Thus, in other embodiments, the invention provides methods 
for identifying a polynucleotide sequence encoding a polypeptide, wherein said 

15 polypeptide may be associated with resistance to development of an infectious disease, 
comprising the steps of: (a) comparing infectious disease-resistant non-human primate 
protein coding sequences to human protein coding sequences, wherein the human protein 
coding sequence is associated with development of the infectious disease; and (b) 
selecting an infectious disease-resistant non-human primate sequence that contains at 

20 least one nucleotide change as compared to the corresponding human sequence, wherein 
the nucleotide change is evolutionarily significant. In another embodiment, the invention 
provides methods for identifying a human polynucleotide sequence encoding a 
polypeptide, wherein said polypeptide may be associated with susceptibility to 
development of an infectious disease, comprising the steps of: (a) comparing human 

25 protein coding sequences to protein-coding polynucleotide sequences of an infectious 
disease-resistant non-human primate, wherein the human protein coding sequence is 
associated with development of the infectious disease; and (b) selecting a human 
polynucleotide sequence that contains at least one nucleotide change as compared to the 
corresponding sequence of an infectious disease-resistant non-human primate, wherein 
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the nucleotide change is evolutionarily significant. 

In the present invention, human sequences to be compared with a homologue from 
an AIDS -resistant non-human primate are selected based on their known or implicated 
association with HIV propagation (Le., replication), dissemination and/or subsequent 

5 progression to AIDS. Such knowledge is obtained, for example, from published literature 
and/or public databases (including sequence databases such as GenBank). Because the 
pathway involved in development of AIDS (including viral replication) involves many 
genes, a number of suitable candidates may be tested using the methods of this invention. 
Table 1 contains a exemplary list of genes to be examined. The sequences are generally 

10 known in the art. 

Table 1: Sample List of Human Genes to be/have been Examined 
Gene Function 



15 eIF-5A 
hPC6A 



initiation factor 
protease 



hPC6B 

20 protease 
P56 lck 

Signal transduction 
FK506-binding protein 

Immunophilin 

25 calnexin 

? 

Bax 

PCD promoter 

bcl-2 

30 apoptosis inhibitor 

lck 

tyrosine kinase 
MAPK (mitogen activated protein kinase) 

protein kinase 

35 CD43 

sialoglycoprotein 

CCR2B 
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chemokine receptor 


CCR3 






chemokine receptor 


Bonzo 






chemokine receptor 


BOB 






chemokine receptor 


GPR1 






chemokine receptor 



10 stromal-derived factor- 1 (SDF-1) 

chemokine 
tumor-necrosis factor-a (TNF- a) 

PCD promoter 

TNF-receptor II (TNFRII) 
15 receptor 
interferon y (IFN- y) 

cytokine 
interleukin 1 a(IL-l a) 

cytokine 

20 interleukin lp(IL-l p) 
cytokine 
interleukin 2 (IL-2) 

cytokine 

interleukin 4 (IL-4) 
25 cytokine 
interleukin 6 (IL-6) 

cytokine 

interleukin 10 (IL- 10) 
cytokine 
30 interleukin 13 (EL-13) 
cytokine 

B7 

signaling protein 
macrophage colony-stimulating factor (M-CSF) 
35 cytokine 

granulocyte-macrophage colony-stimulating factor 

cytokine 

phosphatidylinositol 3-kinase (PI 3-kinase) 

kinase 

40 phosphatidylinositol 4-kinase (PI 4-kinase) 

kinase 

HLA class I a chain 

histocompatibility antigen 
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^2 microglobulin 

lymphocyte antigen 

CD55 

decay-accelerating factor 

5 CD63 

glycoprotein antigen 

CD71 

? 

interferon a (IFN- a) 
10 cytokine 
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CD44 
CD8 



cell adhesion 
glycoprotein 

Genes already examined (13) 



ICAM-1 

20 Immune system 



ICAM-2 
ICAM-3 



Immune system 



Immune system 
25 leukocyte associated function 1 molecule a (LFA-1) 

Immune system 
leukocyte associated function 1 molecule P (LFA-1) 

Immune system 

30 Mac-1 a 

Immune system 
Mac-1 p (equivalent to LFA-1 p) 

Immune system 

DC-SIGN 

35 Immune system 

CD59 

CXCR4 

40 CCR5 

MJP-la 



complement protein 
chemokine receptor 
chemokine receptor 

chemokine 
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MIP-ip 

chemokine 

RANTES 

chemokine 

5 

Aligned protein-coding sequences of human and an AIDS resistant non-human 
primate such as chimpanzee are analyzed to identify nucleotide sequence differences at 
particular sites. The detected sequence changes are generally, and preferably, initially 

10 checked for accuracy as described above. The evolutionarily significant nucleotide 

changes, which are detected by molecular evolution analysis such as the K A /Ks analysis, 
can be further assessed to determine whether the non-human primate gene or the human 
gene has been subjected to positive selection. For example, the identified changes can be 
tested for presence/absence in other AIDS- resistant non-human primate sequences. The 

15 sequences with at least one evolutionarily significant change between human and one 

AIDS-resistant non-human primate can be used as primers for PCR analysis of other non- 
human primate protein-coding sequences, and resulting polynucleotides are sequenced to 
see whether the same change is present in other non-human primates. These comparisons 
allow further discrimination as to whether the adaptive evolutionary changes are unique 

20 to the AIDS-resistant non-human primate (such as chimpanzee) as compared to other 
non-human primates. For example, a nucleotide change that is detected in chimpanzee 
but not other primates more likely represents positive selection on the chimpanzee gene. 
Other non-human primates used for comparison can be selected based on their 
phylogenetic relationships with human. Closely related primates can be those within the 

25 hominoid sublineage, such as chimpanzee, bonobo, gorilla, and orangutan. Non-human 
primates can also be those that are outside the hominoid group and thus not so closely 
related to human, such as the Old World monkeys and New World monkeys. Statistical 
significance of such comparisons maybe determined using established available 
programs, e.g., r-test as used by Messier and Stewart (1997) Nature 385:151-154. 

30 Furthermore, sequences with significant changes can be used as probes in 

genomes from different humans to see whether the sequence changes are shared by more 
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than one individual. For example, certain individuals are slower to progress to AIDS 
("slow progresses") and comparison (a) between a chimpanzee sequence and the 
homologous sequence from the slow-progresser human individual and/or (b) between an 
AEDS-susceptible individual and a slow-progresser individual would be of interest. Gene 
sequences from different human populations can be obtained from databases made 
available by, for example, the human genome diversity project or, alternatively, from 
direct sequencing of PCR-amplified DNA from a number of unrelated, diverse human 
populations. The presence of the identified changes in human slow progresses would 
further indicate the evolutionary significance of the changes. 

As is exemplified herein, the CD59 protein, which has been associated with the 
chimpanzee's resistance to the progression of AIDS, exhibits an evolutionarily significant 
nucleotide change relative to human CD59. CD59 (also known as protectin, lF-5Ag, 
HI 9, HRF20, MACIF, MIRL and P-18) is expressed on peripheral blood leukocytes and 
erythrocytes, and functions to restrict lysis of human cells by complement (Meri et al. 
(1996) Biochem. J. 316:923). More specifically, CD59 acts as an inhibitor of membrane 
attack complexes, which are complement proteins that make hole-like lesions in the cell 
membranes. Thus, CD59 protects the cells of the body from the complement arm of its 
own defense system (Meri et al, supra). The chimpanzee homolog of this protein was 
examined because the human homolog has been implicated in the progression of AIDS in 
infected individuals. It has been shown that CD59 is one of the host cell derived proteins 
that is selectively taken up by HIV virions (Frank et al. (1996) AIDS 10:161 1). 
Additionally, it has been shown that HTV virions that have incorporated host cell CD59 
are protected from the action of complement. Thus, in humans, HIV uses CD59 to 
protect itself from attack by the victim's immune system, and thus to further the course of 
infection. As is theorized in the examples, positively-selected chimpanzee CD59 may 
constitute the adaptive change that inhibits disease progression. The virus may be unable 
to usurp the chimpanzee's CD59 protective role, thereby rendering the vims susceptible to 
the chimpanzee's immune system. 

As is further exemplified herein, the DC-SIGN protein has also been determined 
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to be positively selected in the chimpanzee as compared to humans and gorilla. DC- 
SIGN is expressed on dendritic cells and has been documented to provide a mechanism 
for travel of the HIV-1 virus to the lymph nodes where it infects undifferentiated T cells 
(Geijtenbeek, T.B.H. et al (2000) Cell 100:587-597). Infection of the T cells ultimately 

5 leads to compromise of the immune system and subsequently to full-blown AIDS. The 
HIV-1 virus binds to the extracellular portion of DC-SIGN, and then gains access to the T 
cells via their CD4 proteins. DC-SIGN has as its ligand ICAM-3, which has a very high 
K A /Ks ratio. It may be that the positive selection on chimpanzee ICAM-3 was a result of 
compensatory changes to permit continued binding to DC-SIGN. As is theorized in the 

10 examples, positively-selected chimpanzee DC-SIGN may constitute another adaptive 
change that inhibits disease progression. Upon resolution of the three-dimensional 
structure of chimpanzee DC-SIGN and identification of the mechanism by which HIV-1 
is prevented from binding to DC-SIGN, it may be possible to design drugs to mimic the 
effects of chimpanzee DC-SIGN without disrupting the normal functions of human DC- 

15 SIGN. 

Description of the Breast Enhancement Embodiment fan example of a positively selected 
human trait) 

Relative to non-human primates, female humans exhibit pre-pregnancy, pre- 
20 lactation expanded breast tissue. As is discussed in the Examples, this secondary sex 
characteristic is believed to facilitate evolved behaviors in humans associated with long 
term pair bonds and long-term rearing of infants. One aspect of this invention concerns 
identifying those human genes that have been positively selected in the development of 
enlarged breasts. Specifically, this invention includes a method of determining whether a 
25 human polynucleotide sequence which has been associated with enlarged breasts in 

humans has undergone evolutionarily significant change relative to a non-human primate 
that does not manifest enlarged breasts, comprising: a) comparing the human 
polynucleotide sequence with the corresponding non-human primate polynucleotide 
sequence to identify any nucleotide changes; and b) determining whether the human 
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nucleotide changes are evolutionarily significant. 

It has been found that the human BRCA1 gene, which has been associated with 
normal breast development in humans, has been positively selected relative to the BRCA1 
gene of chimpanzees and other non-human primates. The identified evolutionarily 
5 significant nucleotide changes could be useful in developing agents that can modulate the 
function of the BRCA1 gene or protein. 

Therapeutic compositions that comprise agents 

As described herein, agents can be screened for their capacity to increase or 

10 decrease the effectiveness of the positively selected polynucleotide or polypeptide 

identified according to the subject methods. For example, agents that may be suitable for 
enhancing breast development may include those which interact directly with the BRCA1 
protein or its ligand, or which block inhibitors of BRCA1 protein. Alternatively, an agent 
may enhance breast development by increasing BRCA1 expression. As the mechanism 

15 of BRCA1 is further elucidated, strategies for enhancing its efficacy can be devised. 

In another example, agents that may be suitable for reducing the progression of 
AIDS could include those which directly interact with the human CD59 protein in a 
manner to make the protein unusable to the HIV virion, possibly by either rendering the 
human CDS 9 unsuitable for packing in the virion particle or by changing the orientation 

20 of the protein with respect to the cell membrane (or via some other mechanism). The 
candidate agents can be screened for their capacity to modulate CD59 function using an 
assay in which the agents are contacted with HIV infected cells which express human 
CD59, to determine whether syncytia formation or other indicia of the progression of 
AIDS are reduced. The assay may permit the detection of whether the HIV virion can 

25 effectively pack the CD59 and/or utilize the CD59 to inhibit attack by MAC complexes. 
One agent that may slow AIDS progression is a human CD59 that has been 
modified to have multiple GPI links. As described herein, chimp CD59, which contains 
three GPI links as compared to the single GPI link found in human CD59, slows 
progression of HIV infections in chimps. Preferably, the modified human CD59 contains 
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three GPI links in tandem, 

Another example of an agent that may be suitable for reducing AIDS progression 

is a compound that directly interacts with human DC-SIGN to reduce its capacity to bind 

to HIV-1 and transport it to the lymph nodes. Such an agent could bind directly to the 
5 HIV-1 binding site on DC-SIGN. The candidate agents can be contacted with dendritic 

cells expressing DC-SIGN or with a purified extracellular fragment of DC-SIGN and 

tested for their capacity to inhibit HIV-1 binding. 

Various delivery systems are known in the art that can be used to administer 

agents identified according to the subject methods. Such delivery systems include 
10 aqueous solutions, encapsulation in liposomes, microparticles or microcapsules or 

conjugation to a moiety that facilitates intracellular admission. 

Therapeutic compositions comprising agents may be administered parenterally by 

injection, although other effective administration forms, such as intra- articular injection, 

inhalant mists, orally-active formulations, transdermal iontophoresis or suppositories are 
1 5 also envisioned. The carrier may contain other pharmacologically-acceptable excipients 

for modifying or maintaining the pH, osmolality, viscosity, clarify, color, sterility, 

stability, rate of dissolution, or odor of the formulation. The carrier may also contain 

other pharmacologically-acceptable excipients for modifying or maintaining the stability, 

rate of dissolution, release or absorption of the agent. Such excipients are those 
20 substances usually and customarily employed to formulate dosages for parenteral 

administration in either unit dose or multi-dose form. 

Once the therapeutic composition has been formulated, it may be stored in sterile 

vials as a solution, suspension, gel, emulsion, solid, or dehydrated or lyophilized powder. 

Such formulations may be stored either in a ready to use form or requiring reconstitution 
25 immediately prior to administration. The manner of administering formulations 

containing agents for systemic delivery may be via subcutaneous, intramuscular, 

intravenous, intranasal or vaginal or rectal suppository. Alternatively, the formulations 

may be administered directly to the target organ (e.g., breast). 

The amount of agent which will be effective in the treatment of a particular 



61 



disorder or condition will depend on the nature of the disorder or condition which can be 
determined by standard clinical techniques. In addition, in vitro or in vivo assays may 
optionally be employed to help identify optimal dosage ranges. The precise dose to be 
employed in the formulation will also depend on the route of administration, and the 
5 seriousness or advancement of the disease or condition, and should be decided according 
to the practitioner and each patient's circumstances. Effective doses may be extrapolated 
from dose-response curves derived from in vitro or animal model test systems. For 
example, an effective amount of an agent identified according to the subject methods is 
readily determined by administering graded doses of a bivalent compound of the 

10 invention and observing the desired effect. 

The following examples are provided to further assist those of ordinary skill in the 
art. Such examples are intended to be illustrative and therefore should not be regarded as 
limiting the invention. A number of exemplary modifications and variations are 
described in this application and others will become apparent to those of skill in this art. 

15 Such variations are considered to fall within the scope of the invention as described and 
claimed herein. 



20 EXAMPLES 

EXAMPLE 1: cDNA Library Construction 

A chimpanzee cDNA library is constructed using chimpanzee tissue. Total RNA 
is extracted from the tissue (RNeasy kit, Quiagen; RNAse-free Rapid Total RNA kit, 5 
25 Prime-3 Prime, Inc.) and the integrity and purity of the RNA are determined according to 
conventional molecular cloning methods. Poly A+ RNA is isolated (Mini-Oligo(dT) 
Cellulose Spin Columns, 5 Prime--3 Prime, Inc.) and used as template for the reverse- 
transcription of cDNA with oligo (dT) as a primer. The synthesized cDNA is treated and 
modified for cloning using commercially available kits. Recombinants are then packaged 
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and propagated in a host cell line. Portions of the packaging mixes are amplified and the 
remainder retained prior to amplification. The library can be normalized and the numbers 
of independent recombinants in the library is determined. 

5 EXAMPLE 2: Sequence Comparison 

Suitable primers based on a candidate human gene are prepared and used for PCR 
amplification of chimpanzee cDNA either from a cDNA library or from cDNA prepared 
from mRNA. Selected chimpanzee cDNA clones from the cDNA library are sequenced 
using an automated sequencer, such as an ABI 377. Commonly used primers on the 

1 0 cloning vector such as the M 1 3 Universal and Reverse primers are used to carry out the 
sequencing. For inserts that are not completely sequenced by end sequencing, dye- 
labeled terminators are used to fill in remaining gaps. 

The detected sequence differences are initially checked for accuracy, for example 
by finding the points where there are differences between the chimpanzee and human 

15 sequences; checking the sequence fluorogram (chromatogram) to determine if the bases 
that appear unique to human correspond to strong, clear signals specific for the called 
base; checking the human hits to see if there is more than one human sequence that 
corresponds to a sequence change; and other methods known in the art, as needed. 
Multiple human sequence entries for the same gene that have the same nucleotide at a 

20 position where there is a different chimpanzee nucleotide provides independent support 
that the human sequence is accurate, and that the chimpanzee/human difference is real. 
Such changes are examined using public database information and the genetic code to 
determine whether these DNA sequence changes result in a change in the amino acid 
sequence of the encoded protein. The sequences can also be examined by direct 

25 sequencing of the encoded protein. 

EXAMPLE 3: Molecular Evolution Analysis 

The chimpanzee and human sequences under comparison are subjected toK A /K s 
analysis. In this analysis, publicly available computer programs, such as Li 93 and MA, 
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are used to determine the number of non-synonymous changes per site (K A ) divided by 
the number of synonymous changes per site (Ks) for each sequence under study as 
described above. Full-length coding regions or partial segments of a coding region can be 
used. The higher the K A /Ks ratio, the more likely that a sequence has undergone adaptive 

5 evolution. Statistical significance of K A /Ks values is determined using established 
statistic methods and available programs such as the /-test. 

To further lend support to the significance of a high K A /Ks ratio, the sequence 
under study can be compared in multiple chimpanzee individuals and in other non-human 
primates, e.g., gorilla, orangutan, bonobo. These comparisons allow further 

10 discrimination as to whether the adaptive evolutionary changes are unique to the human 
lineage compared to other non-human primates. The sequences can also be examined by 
direct sequencing of the gene of interest from representatives of several diverse human 
populations to assess to what degree the sequence is conserved in the human species. 

15 EXAMPLE 4: Identification of positively selected ICAM-1, ICAM-2 and ICAM-3 

Using the methods of the invention described herein, the intercellular adhesion 
molecules ICAM-1, ICAM-2 and ICAM-3 have been shown to have been strongly 
positively selected. The ICAM molecules are involved in several immune response 
interactions and are known to play a role in progression to AIDS in HIV infected humans. 

20 The ICAM proteins, members of the Ig superfamily, are ligands for the integrin leukocyte 
associated function 1 molecule (LFA-1). Makgoba et al (1988) Nature 331:86-88. LFA- 
1 is expressed on the surface of most leukocytes, while ICAMs are expressed on the 
surface of both leukocytes and other cell types. Larson et al (1989) /. Cell Biol 
108:703-712. ICAM and LFA-1 proteins are involved in several immune response 

25 interactions, including T-cell function, and targeting of leukocytes to areas of 
inflammation. Larson et al (1989). 

Total RNA was prepared using either the RNeasy® kit (Qiagen), or the RNAse- 
free Rapid Total RNA kit (5 Prime - 3 Prime, Inc.) from primate tissues (chimpanzee 
brain and blood, gorilla blood and spleen, orangutan blood) or from cells harvested from 
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the following B lymphocyte cell lines: CARL (chimpanzee), ROK (gorilla), and PUTI 
(orangutan). mRNA was isolated from total RNA using the Mini-Oligo(dT) Cellulose 
Spin Columns (5 Prime - 3 Prime, Inc.). cDNA was synthesized from mRNA with oligo 
dT and/or random priming using the cDNA Synthesis Kit (Stratagene®). The protein- 
coding region of the primate ICAM-1 gene was amplified from cDNA using primers 
(concentration^ 00 nmole/ul) designed by hand from the published human sequence. 
PCR conditions for ICAM-1 amplification were 94°C initial pre-melt (4 min), followed 
by 35 cycles of 94°C (15 sec), 58°C (1 min 15 sec), 72°C (1 min 15 sec), and a final 72°C 
extension for 10 minutes. PCR was accomplished using Ready-to-Go™ PCR beads 
(Amersham Pharmacia Biotech) in a 50 microliter total reaction volume. Appropriately- 
sized products were purified from agarose gels using the QiaQuick® Gel Extraction kit 
(Qiagen). Both strands of the amplification products were sequenced directly using the 
Big Dye Cycle Sequencing Kit and analyzed on a 373 A DNA sequencer (ABI 
BioSystems). 

Comparison of the protein-coding portions of the human, gorilla (Gorilla gorilla), 
and orangutan (Pongo pygmaeus) ICAM-1 genes to that of the chimpanzee yielded 
statistically significant K A /K S ratios (Table 2). The protein-coding portions of the human 
and chimpanzee ICAM-1 genes were previously published and the protein-coding 
portions of gorilla (Gorilla gorilla), and orangutan (Pongo pygmaeus) ICAM-1 genes are 
shown in Figures 3 and 4, respectively. 

For this experiment, pairwise K A /K S ratios were calculated for the mature protein 
using the algorithm of Li (1985; 1993). Statistically significant comparisons (determined 
by Mests) are shown in bold. Although the comparison to gorilla and human was 
sufficient to demonstrate that chimpanzee ICAM-1 has been positively-selected, the 
orangutan ICAM-1 was compared as well, since the postulated historical range of gorillas 
in Africa suggests that gorillas could have been exposed to the HIV-1 virus. Nowak and 
Paradiso (1983) Walker's Mammals of the World (Baltimore, MD, The Johns Hopkins 
University Press). The orangutan, however, has always been confined to Southeast Asia 
and is thus unlikely to have been exposed to HIV over an evolutionary time frame. 
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(Nowak and Paradiso, 1983) (Gorillas are most closely-related to humans and 
chimpanzees, while orangutans are more distantly-related.) 



Table 2. K A /Ks Ratios: ICAM-1 Whole Protein Comparisons 



Species Compared 


K A /Ks Ratio 


Chimpanzee to Human 


2.1 (P < 0.01) 


Chimpanzee to Gorilla 


1.9 (P < 0.05) 


Chimpanzee to Orangutan 


1.4 (P < 0.05) 


Human to Gorilla 


1.0 


Human to Orangutan 


0.87 


Gorilla to Orangutan 


0.95 



5 

Even among those proteins for which positive selection has been demonstrated, 



few show K A /K S ratios as high as these ICAM-1 comparisons. Lee and Vacquier (1992) 
Biol. Bull. 182:97-104; Swanson and Vacquier (1995) Proc. Natl. Acad. Sci. USA 
92:4957-4961; Messier and Stewart (1997); Sharp (1997) Nature 385:111-112. The 

10 results are consistent with strong selective pressure resulting in adaptive changes in the 
chimpanzee ICAM-1 molecule. 

The domains (Dl and D2) of the ICAM-1 molecule which bind to LFA-1 have 
been documented. Staunton et al. (1990). Cell 61:243-254. Pairwise K A /K S comparisons 
between primate ICAM-1 genes. K A /K S ratios were calculated for domains Dl and D2 

15 only, using the algorithm of Li (1985; 1993) (Table 3). Statistically significant 
comparisons (determined by f-tests) are shown in bold. The very high, statistically 
significant K A /K S ratios for domains Dl and D2 suggest that these regions of the protein 
were very strongly positively-selected. These regions of chimpanzee ICAM-1 display 
even more striking K A /K S ratios (Table 3) than are seen for the whole protein 

20 comparisons, thus suggesting that the ICAM-l/LFA-1 interaction has been subjected to 
unusually strong selective pressures. 

Table 3. K A /Ks Ratios: Domains D1+D2 of ICAM-1 
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Species Compared 


K A /Ks Ratio 


Chimpanzee to Human 


3.1 (P<0.01) 


Chimpanzee to Gorilla 


2.5 (P < 0.05) 


Chimpanzee to Orangutan 


1.5 (P < 0.05) 


Human to Gorilla 


1.0 


Human to Orangutan 


0.90 


Gorilla to Orangutan 


1.0 



EXAMPLE 5: Characterization of ICAM-L ICAM-2 and ICAM-3 positively 
selected sequences 

A sequence identified by the methods of this invention may be further tested and 
5 characterized by cell transfection experiments. For example, human cells in culture, 
when transfected with a chimpanzee polynucleotide identified by the methods described 
herein (such as ICAM-1 (or ICAM-2 or ICAM-3); see below), could be tested for reduced 
viral dissemination and/or propagation using standard assays in the art, and compared to 
control cells. Other indicia may also be measured, depending on the perceived or 

10 apparent functional nature of the polynucleotide/polypeptide to be tested. For example, in 
the case of ICAM-1 (or ICAM-2 or ICAM-3), syncytia formation maybe measured and 
compared to control (untransfected) cells. This would test whether the resistance arises 
from prevention of syncytia formation in infected cells. 

Cells which are useful in characterizing sequences identified by the methods of 

1 5 this invention and their effects on cell-to-cell infection by HIV-1 are human T-cell lines 
which are permissive for infection with HIV-1, including, e.g., H9 and HUT78 cell lines, 
which are available from the ATCC 

For cell transfection assays, ICAM-1 (or ICAM-2 or ICAM-3) cDNA (or any 
cDNA identified by the methods described herein) can be cloned into an appropriate 

20 expression vector. To obtain maximal expression, the cloned ICAM-1 (or ICAM-2 or 

ICAM-3) coding region is operably linked to a promoter which is active in human T cells, 
such as, for example, an IL-2 promoter. Alternatively, an ICAM-1 (or ICAM-2 or ICAM- 
3) cDNA can be placed under transcriptional control of a strong constitutive promoter, or 
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an inducible promoter. Expression systems are well known in the art, as are methods for 
introducing an expression vector into cells. For example, an expression vector 
comprising an ICAM-1 (or ICAM-2 or ICAM-3) cDNA can be introduced into cells by 
DEAE-dextran or by electroporation, or any other known method. The cloned ICAM-1 

5 (or ICAM-2 or ICAM-3) molecule is then expressed on the surface of the cell 

Determination of whether an ICAM-1 (or ICAM-2 or ICAM-3) cDNA is expressed on 
the cell surface can be accomplished using antibody(ies) specific for ICAM-1 (or ICAM-2 
or ICAM-3). In the case of chimpanzee ICAM-1 (or ICAM-2 or ICAM-3) expressed on 
the surface of human T cells, an antibody which distinguishes between chimpanzee and 

10 human ICAM-1 (or ICAM-2 or ICAM-3) can be used. This antibody can be labeled with 
a detectable label, such as a fluorescent dye. Cells expressing chimpanzee ICAM-1 (or 
ICAM-2 or ICAM-3) on their surfaces can be detected using fluorescence-activated cell 
sorting and the anti-ICAM-1 (or ICAM-2 or ICAM-3) antibody appropriately labeled, 
using w r ell-established techniques. 

1 5 Transfected human cells expressing chimpanzee IC AM- 1 (or ICAM-2 or ICAM- 

3) on their cell surface can then be tested for syncytia formation, and/or for HIV 
replication, and/or for number of cells infected as an index of cell-to-cell infectivity. The 
chimpanzee ICAM-1 (or ICAM-2 or ICAM-3)-expressing cells can be infected with HIV- 
1 at an appropriate dose, for example tissue culture infectious dose 50, i.e., a dose which 

20 can infect 50% of the cells. Cells can be plated at a density of about 5 x 10 5 cells/ml in 
appropriate tissue culture medium, and, after infection, monitored for syncytia formation, 
and/or viral replication, and/or number of infected cells in comparison to control, 
uninfected cells. Cells which have not been transfected with chimpanzee ICAM-1 (or 
ICAM-2 or ICAM-3) also serve as controls. Syncytia formation is generally observed in 

25 HIV- 1 -infected cells (which are not expressing chimpanzee ICAM-1 (or ICAM-2 or 
ICAM-3)) approximately 10 days post-infection. 

To monitor HIV replication, cell supernatants can be assayed for the presence and 
amount of p24 antigen. Any assay method to detect p24 can be used, including, for 
example, an ELISA assay in which rabbit anti-p24 antibodies are used as capture 
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antibody, biotinylated rabbit anti-p24 antibodies serve as detection antibody, and the 
assay is developed with avidin-horse radish peroxidase. To determine the number of 
infected cells, any known method, including indirect immunofluorescence methods, can 
be used. In indirect immunofluorescence methods, human HIV-positive serum can be 
5 used as a source of anti-HIV antibodies to bind to infected cells. The bound antibodies 
can be detected using FITC-conjugated anti-human IgG, the cells visualized by 
fluorescence microscopy and counted. 

Another method for assessing the role of a molecule such as ICAM-1 (or ICAM-2 
or ICAM-3) involves successive infection of cells with HIV. Human cell lines, preferably 
10 those that do not express endogenous ICAM (although cell lines that do express 

endogenous ICAM may also be used), are transfected with either human or chimpanzee 
ICAM -1 or -2 or -3. In one set of experiments, HIV is collected from the supernatant of 
HIV-infected human ICAM-1 (or ICAM-2 or ICAM-3)-expressing cells and used to 
infect chimpanzee ICAM-1 (or ICAM-2 or ICAM-3)-expressing cells or human ICAM-1 
1 5 (or ICAM-2 or ICAM-3)-expressing cells. Initial infectivity, measured as described 

above, of both the chimpanzee ICAM-1 (or ICAM-2 or ICAM-3)- and the human ICAM- 
1 (or ICAM-2 or ICAM-3)-expressing cells would be expected to be high. After several 
rounds of replication, cell to cell infectivity would be expected to decrease in the 
chimpanzee ICAM-1 (or ICAM-2 or ICAM-3) expressing cells, if chimpanzee ICAM-1 
20 (or ICAM-2 or ICAM-3) confers resistance. In a second set of experiments, HIV is 
collected from the supernatant of HIV-infected chimpanzee ICAM-1 (or ICAM-2 or 
ICAM-3)-expressing cells, and used to infect human ICAM-1 (or ICAM-2 or ICAM-3)- 
expressing cells. In this case, the initial infectivity would be expected to be much lower 
than in the first set of experiments, if ICAM-1 (or ICAM-2 or ICAM-3) is involved in 
25 susceptibility to HIV progression. After several rounds of replication, the cell to cell 
infectivity would be expected to increase. 

The identified human sequences can be used in establishing a database of 
candidate human genes that may be involved in conferring, or contributing to, AIDS 
susceptibility or resistance. Moreover, the database not only provides an ordered 
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collection of candidate genes, it also provides the precise molecular sequence differences 
that exist between human and an AIDS-resistant non-human primate (such as 
chimpanzee) and thus defines the changes that underlie the functional differences. 

5 EXAMPLE 6: Molecular Modeling of ICAM-1 and ICAM-3 

Modeling of the three-dimensional structure of ICAM-1 and ICAM-3 has 
provided additional evidence for the role of these proteins in explaining chimpanzee 
resistance to AIDS progression. 

In the case of ICAM-1, 5 of the 6 amino acid replacements that are unique to the 
1 0 chimpanzee lineage are immediately adjacent (i.e., physically touching) to those amino 
acids identified by mutagenic studies as critical to LFA-1 binding. These five amino acid 
replacements are human LI 8 to chimp Q18, human K29 to chimp D29, human P45 to 
chimp G45, human R49 to chimp W49, and human E171 to chimp Q171. This 
positioning cannot be predicted from the primary structure (i.e., the actual sequence of 
1 5 amino acids). None of the amino acid residues critical for binding has changed in the 
chimpanzee ICAM-1 protein. 

Such positioning argues strongly that the chimpanzee ICAM-1 protein's basic 
function is unchanged between humans and chimpanzees; however, evolution has 
wrought fine-tuned changes that may help confer upon chimpanzees their resistance to 
20 progression of AIDS. The nature of the amino acid replacements is being examined to 
allow exploitation of the three-dimensional structural information for developing agents 
for therapeutic intervention. Strikingly, 4 of the 5 chimpanzee residues are adjacent to 
critical binding residues that have been identified as N-linked glycosylation sites. This 
suggests that differences exist in binding constants (to LFA-1) for human and chimpanzee 
25 ICAM- 1 . These binding constants are being determined. Should the binding constants 
prove lower in chimpanzee ICAM-1, it is possible to devise small molecule agents to 
mimic (by way of steric hindrance) the change in binding constants as a potential 
therapeutic strategy for HIV-infected humans. Similarly, stronger binding constants, if 
observed for chimpanzee ICAM-1, will suggest alternative strategies for developing 
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therapeutic interventions for HIV-1 infected humans. 

In the case of ICAM-3, a critical amino acid residue replacement from proline 
(observed in seven humans) to glutamine (observed in three chimpanzees) is predicted 
from our modeling studies to significantly change the positional angle between domains 2 
5 and 3 of human and chimpanzee ICAM-3 . The human protein displays an acute angle at 
this juncture. Klickstein, et ai, 1996 J. Biol. Chem. 27:239 20-27. Loss of this sharp 
angle (bend) is predicted to render chimpanzee ICAM-3 less easily packaged into HIV-1 
virions (In infected humans, after ICAMs are packaged into HIV virions, cell-to-cell 
infectivity dramatically increases. Barbeau, B. et al, 1998 J. Virol. 72:7125-7136). This 
1 0 failure to easily package chimp ICAM-3 into HIV virions could then prevent the increase 
in cell-to-cell infectivity seen in infected humans. This would then account for 
chimpanzee resistance to AIDS progression. 

A small molecule therapeutic intervention whereby binding of a suitably-designed 
small molecule to the human proline residue causes (as a result of steric hindrance) the 
15 human ICAM-1 protein to mimic the larger (i.e., less-acute) angle of chimpanzee ICAM- 
3 is possible. Conservation between the 2 proteins of the critical binding residues (and 
the general resemblance of immune responses between humans and chimpanzees) argues 
that alteration of this angle will not compromise the basic function of human ICAM-3. 
However, the human ICAM-3 protein would be rendered resistant to packaging into HTV 
20 virions, thus mimicking (in HIV-1 infected humans) the postulated pathway by which 
infected chimpanzees resist progression to AIDS. 

Essentially the same procedures were used to identify positively selected 
chimpanzee ICAM-2 and ICAM-3 (see Table 4). The ligand binding domain of ICAM-1 
has been localized as exhibiting especially striking positive selection in contrast to 
25 ICAMs -2 and -3, for which positive selection resulted in amino acid replacements 
throughout the protein. Thus, this comparative genomic analysis reveals that positive 
selection on ICAMs in chimpanzees has altered the proteins' primary structure, for 
example, in important binding domains. These alterations may have conferred resistance 
to AIDS progression in chimpanzees. 
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Table 4. K A /Ks Ratios: ICAM-2 and 3 Whole Protein Comparisons 



Species Compared 


Ka/Ks Ratio 


Chimpanzee to Human ICAM-2 


2.1 (P<0.01) 


Chimpanzee to Human ICAM-3 


3.7 (P < 0.01) 



Binding of ICAM-1, -2, and -3 has been demonstrated to play an essential role in 
the formation of syncytia (i.e., giant, multi-nucleated cells) in HIV-infected cells in vitro. 
5 Pantaleo et al. (1991) J. Ex. Med. 173:51 1-514. Syncytia formation is followed by the 
depletion of CD + cells in vitro. Pantaleo et al (1991); Levy (1993) Microbiol. Rev. 
57:183-189; Butini et al (1994) Eur. J. Immunol. 24:2191-2195; Finkel and Banda 
(1994) Curr. Opin. Immunol. 6:605-615. Although syncytia formation is difficult to 
detect in vivo, clusters of infected cells are seen in lymph nodes of infected individuals. 
10 Pantaleo et al, (1993) N. Eng. J. Med. 328:327-335; Finkel and Banda (1994); 

Embretson et al. (1993) Nature 362:359-362; Pantaleo et al. (1993) Nature 362:355-358. 
Syncytia may simply be scavenged from the body too quickly to be detected. Fouchier et 
al (1996) Virology 219:87-95. Syncytia-mediated loss of CD4 + cells in vivo has been 
speculated to occur; this could contribute directly to compromise of the immune system, 
1 5 leading to opportunistic infection and full-blown AIDS. Sodrosky et al. (1 986) Nature 
322:470-474; Hildreth and Orentas (1989) Science 244:1075-1078; Finkel and Banda 
(1994). Thus critical changes in chimpanzee ICAM-1, ICAM-2 or ICAM-3 may deter 
syncytia formation in chimpanzee and help explain chimpanzee resistance to AIDS 
progression. Because of the polyfunctional nature of ICAMs, these positively selected 
20 changes in the ICAM genes may additionally confer resistance to other infectious diseases 
or may play a role in other inflammatory processes that may also be of value in the 
development of human therapeutics. The polypeptide sequence alignments of ICAM-1, - 
2, and -3 are shown in Figures 5, 6, and 7, respectively. 

25 EXAMPLE 7: Identifying Positive Se lection of MIP-1 a 

MlP-la is a chemokine that has been shown to suppress HTV-1 replication in 
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human cells in vitro (Cocchi, F. etaL, 1995 Science 270:1811-1815). The chimpanzee 
homologue of the human MlP-la gene was PCR-amplified and sequenced. Calculation 
of the K A /Ks ratio (2.1, P<0.05) and comparison to the gorilla homologue reveals that the 
chimpanzee gene has been positively-selected. As for the other genes discussed herein, 
5 the nature of the chimpanzee amino acid replacements is being examined to determine 
how to exploit the chimpanzee protein for therapeutic intervention. 

EXAMPLE 8: Identifying Positive Selection of 17-3-Hydroxysteroid Dehydrogenase 

Using the methods of the present invention, a chimpanzee gene expressed in brain 
has been positively-selected (K A /Ks =1.6) as compared to its human homologue 
10 (GenBank Acc. # X87176) has been identified. The human gene, 17-P hydroxysteroid 
*f dehydrogenase type IV, codes for a protein known to degrade the two most potent 

a estrogens, P-estradiol, and 5-diol (Adamski, J. et al. 1995 Biochem J. 31 1:437-443). 

fy Estrogen-related cancers (including, for example, breast and prostate cancers) account for 

! 2 some 40% of human cancers. Interestingly, reports in the literature suggest that 

15 chimpanzees are resistant to tumorigenesis, especially those that are estrogen-related, 
p This protein may have been positively-selected in chimpanzees to allow more efficient 

degradation of estrogens, thus conferring upon chimpanzees resistance to such cancers. If 
y so, the specific amino acid replacements observed in the chimpanzee protein may supply 

M important information for therapeutic intervention in human cancers. 

20 

EXAMPLE 9: cDNA Library construction for Chimpanzee Brain Tissue 

A chimpanzee brain cDNA library is constructed using chimpanzee brain tissue. 
The chimpanzee brain tissue can be obtained after natural death so that no killing of an 
animal is necessary for this study. In order to increase the chance of obtaining intact 
25 mRNAs expressed in brain, however, the brain is obtained as soon as possible after the 
animal's death. Preferably, the weight and age of the animal are determined prior to 
death. The brain tissue used for constructing a cDNA library is preferably the whole 
brain in order to maximize the inclusion of mRNA expressed in the entire brain. Brain 
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tissue is dissected from the animal following standard surgical procedures. 

Total RNA is extracted from the brain tissue and the integrity and purity of the 
RNA are determined according to conventional molecular cloning methods. Poly A+ 
RNA is selected and used as template for the reverse-transcription of cDNA with oligo 
5 (dT) as a primer. The synthesized cDNA is treated and modified for cloning using 
commercially available kits. Recombinants are then packaged and propagated in a host 
cell line. Portions of the packaging mixes are amplified and the remainder retained prior 
to amplification. The library can be normalized and the numbers of independent 
recombinants in the library is determined. 

10 

EXAMPLE 10; Sequence Comparison of Chimpanzee and H uman Brain cDNA 

Randomly selected chimpanzee brain cDNA clones from the cDNA library are 
sequenced using an automated sequencer, such as the ABI 377. Commonly used primers 
on the cloning vector such as the M13 Universal and Reverse primers are used to carry 

1 5 out the sequencing. For inserts that are not completely sequenced by end sequencing, 
dye-labeled terminators are used to fill in remaining gaps. 

The resulting chimpanzee sequences are compared to human sequences via 
database searches, e.g., BLAST searches. The high scoring "hits," i.e., sequences that 
show a significant (e.g., >80%) similarity after BLAST analysis, are retrieved and 

20 analyzed. The two homologous sequences are then aligned using the alignment program 
CLUSTAL V developed by Higgins et al. Any sequence divergence, including 
nucleotide substitution, insertion and deletion, can be detected and recorded by the 
alignment. 

The detected sequence differences are initially checked for accuracy by finding the 
25 points where there are differences between the chimpanzee and human sequences; 

checking the sequence fluorogram (chromatogram) to determine if the bases that appear 
unique to human correspond to strong, clear signals specific for the called base; checking 
the human hits to see if there is more than one human sequence that corresponds to a 
sequence change; and other methods known in the art as needed. Multiple human 
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sequence entries for the same gene that have the same nucleotide at a position where there 
is a different chimpanzee nucleotide provides independent support that the human 
sequence is accurate, and that the chimpanzee/human difference is real. Such changes are 
examined using public database information and the genetic code to determine whether 
5 these DNA sequence changes result in a change in the amino acid sequence of the 
encoded protein. The sequences can also be examined by direct sequencing of the 
encoded protein. 

EXAMPLE 11: Molecular Evolution Analysis of Human Brain Sequences Relative 

10 to Other Primates 

The chimpanzee and human sequences under comparison are subjected to K A /K S 
analysis. In this analysis, publicly available computer programs, such as Li 93 and INA, 
are used to determine the number of non-synonymous changes per site (K A ) divided by 
the number of synonymous changes per site (K s ) for each sequence under study as 
1 5 described above. This ratio, K A /Ks, has been shown to be a reflection of the degree to 
which adaptive evolution, i.e., positive selection, has been at work in the sequence under 
study. Typically, full-length coding regions have been used in these comparative 
analyses. However, partial segments of a coding region can also be used effectively. The 
higher the K A /K S ratio, the more likely that a sequence has undergone adaptive evolution. 
20 Statistical significance of K A /Ks values is determined using established statistic methods 
and available programs such as the Mest. Those genes showing statistically high K A /K S 
ratios between chimpanzee and human genes are very likely to have undergone adaptive 
evolution. 

To further lend support to the significance of a high K A /Ks ratio, the sequence 
25 under study can be compared in other non-human primates, e.g., gorilla, orangutan, 
bonobo. These comparisons allow further discrimination as to whether the adaptive 
evolutionary changes are unique to the human lineage compared to other non-human 
primates. The sequences can also be examined by direct sequencing of the gene of 
interest from representatives of several diverse human populations to assess to what 
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degree the sequence is conserved in the human species. 

EXAMPLE 12: Further Sequence Characterization of Selected Human Brain 
Sequences 

5 Human brain nucleotide sequences containing evolutionarily significant changes 

are further characterized in terms of their molecular and genetic properties, as well as 
their biological functions. The identified coding sequences are used as probes to perform 
in situ mRNA hybridization that reveals the expression pattern of the gene, either or both 
in terms of what tissues and cell types in which the sequences are expressed, and when 

1 0 they are expressed during the course of development or during the cell cycle. Sequences 
that are expressed in brain may be better candidates as being associated with important 
human brain functions. Moreover, the putative gene with the identified sequences are 
subjected to homologue searching in order to determine what functional classes the 
sequences belong to. 

15 Furthermore, for some proteins, the identified human sequence changes may be 

useful in estimating the functional consequence of the change. By using such criteria a 
database of candidate genes can be generated. Candidates are ranked as to the likelihood 
that the gene is responsible for the unique or enhanced abilities found in the human brain 
compared to chimpanzee or other non-human primates, such as high capacity information 

20 processing, storage and retrieval capabilities, language abilities, as well as others. In this 
way, this approach provides a new strategy by which such genes can be identified. Lastly, 
the database not only provides an ordered collection of candidate genes, it also provides 
the precise molecular sequence differences that exist between human and chimpanzee 
(and other non-human primates), and thus defines the changes that underlie the functional 

25 differences. 

In some cases functional differences are evaluated in suitable model systems, 
including, but not limited to, in vitro analysis such as indicia of long term potentiation 
(LTP), and use of transgenic animals or other suitable model systems. These will be 
immediately apparent to those skilled in the art. 
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EXAMPLE 13: Identification of Positive Selection in a H uman Tyrosine Kinase 
Gene 

Using the methods of the present invention, a human gene (GenBank Acc.# 
5 ABO 14541), expressed in brain has been identified, that has been positively-selected as 
compared to its gorilla homologue. This gene, which codes for a tyrosine kinase, is 
homologous to a well-characterized mouse gene (GenBank Acc.# AF011908) whose gene 
product, called AATYK, is known to trigger apoptosis (Gaozza, E. et al. 1997 Oncogene 
15:3127-3135). The literature suggests that this protein controls apoptosis in the 
1 0 developing mouse brain (thus, in effect, "sculpting" the developing brain). The AATYK- 
induced apoptosis that occurs during brain development has been demonstrated to be 
necessary for normal brain development. 

There is increasing evidence that inappropriate apoptosis contributes to the 
pathology of human neurodegenerative diseases, including retinal degeneration, 
1 5 Huntington's disease, Alzheimer's disease, Parkinson' s disease and spinal muscular 

atrophy, an inherited childhood motoneuron disease. On the other hand in neural tumour 
cells, such as neuroblastoma and medulloblastoma cells, apoptotic pathways may be 
disabled and the cells become resistant to chemotherapeutic drugs that kill cancer cells by 
inducing apoptosis. A further understanding of apoptosis pathways and the function of 
20 apoptosis genes should lead to a better understanding of these conditions and permit the 
use of AATYKI in diagnosis of such conditions. 

Positively-selected human and chimpanzee AATYK may constitute another 
adaptive change that has implications for disease progression. Upon resolution of the 
three-dimensional structure of human and chimpanzee AATYK, it may be possible to 
25 design drugs to modulate the function of AATYK in a desired manner without disrupting 
any of the normal functions of human AATTK. 

It has been demonstrated that mouse AATYK is an active, non-receptor, cytosolic 
kinase which induces neuronal differentiation in human adrenergic neuroblastoma 
(NB):SH-SY5Y cells. AATYK also promotes differentiation induced by other agents, 
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including all-trans retinoic acid (RA), 12-O-Tetradecanoyl phorbol 13-acetate (TP A) and 
IGF-L Raghunath, et ah, Brain Res Mol Brain Res. (2000) 77: 1 5 1-62. In experiments 
with rats, it was found that the AATYK protein was expressed in virtually all regions of 
the adult rat brain in which neurons are present, including olfactory bulb, forebrain, 

5 cortex, midbrain, cerebellum and pons. Immunohistochemical labeling of adult brain 
sections showed the highest levels of AATYK expression in the cerebellum and olfactory 
bulb. Expression of AATYK was also up-regulated as a function of retinoic acid-induced 
neuronal differentiation of pl9 embryonal carcinoma cells, supporting a role for this 
protein in mature neurons and neuronal differentiation. Baker, et al., Oncogene (2001) 

10 20:1015-21. 

Nicolini, et al., Anticancer Res (1998) 18:2477-81 showed that retinoic acid (RA) 
differentiated SH-SY5Y cells were a suitable and reliable model to test the neurotoxicity 
of chemotherapic drugs without the confusing effects of the neurotrophic factors 
commonly used to induce neuronal differentiation. The neurotoxic effect and the course 

15 of the changes is similar to that observed in clinical practice and in in vivo experimental 
models. Thus, the model is proposed as a screening method to test the neurotoxicity of 
chemotherapy drugs and the possible effect of neuroprotectant molecules and drugs. 
Similarly, AATYK differentiated SYSY-5Y cells could be used as a model for screening 
chemotherapeutic drugs and possible side effects of neuroprotectant molecules and drugs. 

20 It has also been shown that AATYK mRNA is expressed in neurons throughout 

the adult mouse brain. AATYK possessed tyrosine kinase activity and was 
autophosphorylated when expressed in 293 cells. AATYK mRNA expression was 
rapidly induced in cultured mouse cerebellar granule cells during apoptosis induced by 
KC1. The number of apoptotic granule cells overexpressing wild-type AATYK protein 

25 was significantly greater than the number of apoptotic granule cells overexpressing a 

mutant AATYK that lacked tyrosine kinase activity. These findings suggest that through 
its tyrosine kinase activity, AATYK is also involved in the apoptosis of mature neurons. 
Tomomura, et al, Oncogene (2001) 20(9): 1022-32. 

The tyrosine kinase domain of AATYK protein is highly conserved between 
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mouse, chimpanzee, and human (as are most tyrosine kinases). Interestingly, however, 
the region of the protein to which signaling proteins bind has been positively-selected in 
humans, but strongly conserved in both chimpanzees and mice. The region of the human 
protein to which signaling proteins bind has not only been positively-selected as a result 

5 of point nucleotide mutations, but additionally displays duplication of several src 
homology 2 (SH2) binding domains that exist only as single copies in mouse and 
chimpanzee. This suggests that a different set of signaling proteins may bind to the 
human protein, which could then trigger different pathways for apoptosis in the 
developing human brain compared to those in mice and chimpanzees. Such a gene thus 

10 may contribute to unique or enhanced human cognitive abilities. Human AATYK has 
been mapped on 25.3 region of chromosome 17. Seki, et al., J Hum Genet (1999) 
44:141-2. 

Chimpanzee DNA was sequenced as part of a high-throughput sequencing project 
on a MegaBACE 1000 sequencer (AP Biotech). DNA sequences were used as query 

15 sequences in a BLAST search of the GenBank database. Two random chimpanzee 
sequences, termed stch856 and stch610, returned results for two genes in the non- 
redundant database of GenBank: NM_004920 (human apoptosis-associated tyrosine 
kinase, AATYK) and ABO 14541 (human KIAA641, identical nucleotide sequence to 
NM__004920), shown in Figure 14 A, and also showed a high K A /Ks ratio compared to 

20 these human sequences. Primers were designed for PCR and sequencing of AATYK. 
Sequence was obtained for the 3 prime end of this gene in chimp and gorilla. The 5 
prime end of the gene was difficult to amplify, and no sequence was confirmed in human 
and gorilla. The human AATYK gene (SEQ ID NO: 14) has a coding region of 3624 bp 
(nucleotides 413-4036 of SEQ ID NO:14), and codes for a protein of 1207 amino acids 

25 (SEQ ID NO: 16). 1 809 bp were sequenced in both chimp and gorilla. See Figure 1 5 A 
and 15B. The partial sequences (SEQ ID NO: 17 and SEQ ID NO: 18) did not include the 
start or stop codons, although they were very close to the stop codon on the 3 prime end 
(21 codons away). These sequences correspond to nucleotides 2170-3976 or 2179-3988 
of the corresponding human sequences taking into account the gaps described below. 
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There were also several pairs of amino acid insertions/deletions among chimp, 
human and gorilla in the coding region. The following sequences are in reading frame: 



15 



Chimp 




GGTGAGGGCCCCGGCCCCGGGCCC 




(SEQ 


ID 


NO: 


19 


Human 


2819 


GGTGAGGGC :::::: CCCGGGCCC 


2836 


(SEO 


ID 


NO : 


20 


Gorilla 




GGCGAGGGC :::::: CCCGGGCCC 




(SEQ 


ID 


NO: 


21 


Chimp 




CTGGAGGCTGAGGCCGAGGCCGAG 




(SEQ 


ID 


NO: 


22 


Human 


2912 


CTCGAGGCT :::::: GAGGCCGAG 


2929 


(SEQ 


ID 


NO: 


23 


Gorilla 




CTGGAGGCT :::::: GAGGCCGAG 




(SEQ 


ID 


NO: 


24 


Chimp 




CCCACGCCC :::::: GCTCCCTTC 




(SEQ 


ID 


NO: 


25 


Human 


3890 


CCCACGCCCACGCCCGCTCCCTTC 


3913 


(SEQ 


ID 


NO: 


26 


Gorilla 




CCCACGCCC: : : : : : GCTCCCTTC 




(SEQ 


ID 


NO: 


27 


Chimp 




CCCACGTCCACGTCCCGCTTCTCC 




(SEQ 


ID 


NO: 


28 


Human 


3938 


CCCACGTCC :::::: CGCTTCTCC 


3955 


(SEQ 


ID 


NO: 


29 


Gorilla 




CCCACGTCC :::::: CGCTTCTCC 




(SEQ 


ID 


NO: 


30 



20 Each of these insertions/deletions affected two amino acids and did not change the 

reading frame of the sequence. Sliding window K A /K S for chimp to human, chimp to 
gorilla, and human to gorilla, excluding the insertion/deletion regions noted above, 
showed a high Ka/Ks ratio for some areas. See Table 9. 

The highest Ka/Ks ratios are human to gorilla and chimp to gorilla, suggesting 

25 that both the human and chimp gene have undergone selection, and is consistent with the 
idea that the two species share some enhanced cognitive abilities relative to the other 
great apes (gorillas, for example). Such data bolsters the view that this gene may play a 
role with regard to enhanced cognitive functions. It should also be noted that in general, 
the human-containing pairwise comparisons are higher than the analogous chimp- 

30 containing comparisons. 
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EXAMPLE 14: Positively Selected Human BRCA1 Gene 

Comparative evolutionary analysis of the BRCA1 genes of several primate species 
has revealed that the human BRCA1 gene has been subjected to positive selection. 
5 Initially, 1141 codons of exon 1 1 of the human and chimpanzee BRCA1 genes (Hacia et 
al (1998) Nature Genetics 18:155-158) were compared and a strikingly high K A /Ks ratio, 
3.6, was found when calculated by the method of Li (Li (1993) J. Mol Evol 36:96-99; Li 
et al (1985) Mol Biol Evol 2:150-174). In fact, statistically significant elevated ratios 
were obtained for this comparison regardless of the particular algorithm used (see Table 

10 5 A). Few genes (or portions of genes) have been documented to display ratios of this 

magnitude (Messier et al (1997) Nature 385:151-154; Endo et al (1996) Mol Biol Evol 
13:685-690; and Sharp (1997) Nature 385:1 1 1-1 12). We thus chose to sequence the 
complete protein-coding region (5589 bp) of the chimpanzee BRCA1 gene, in order to 
compare it to the full-length protein-coding sequence of the human gene. In many cases, 

1 5 even when positive selection can be shown to have operated on limited regions of a 

particular gene, K A /K S analysis of the full-length protein-coding sequence fails to reveal 
evidence of positive selection (Messier et al (1997), supra). This is presumably because 
the signal of positive selection can be masked by noise when only small regions of a gene 
have been positively selected, unless selective pressures are especially strong. However, 

20 comparison of the full-length human and chimpanzee BRCA1 sequences still yielded 

Ka/Ks ratios in excess of one, by all algorithms we employed (Table 5 A). This suggests 
that the selective pressure on BRCA1 was intense. A sliding-window K A /Ks analysis was 
also performed, in which intervals of varying lengths (from 150 to 600 bp) were 
examined, in order to determine the pattern of selection within the human BRCA1 gene. 

25 This analysis suggests that positive selection seems to have been concentrated in exon 1 1 . 



82 



Table 5A: Human-Chimpanzee K A /Ks Comparisons 



Method 


K A /Ks (exon 11) 


K A /Ks (full-length) 


Li (1993) J. Mol Evol 36:96; Li 
etal (1985) Mol Biol Evol 
2:150 


3.6*** 


2.3* 


Ina Y.{\99S)1 Mol Evol 
40:190 


3 3** 


2.1* 


Kumar et al, MEGA: Mol Evol 
Gen. Anal (PA St. Univ, 1993) 


2.2* 


1.2 



Table 5B: K A /Ks for Exon 11 of BRCA1 from Additional Primates 



Comparison 


K A 


Ks 


K A /K S 



Human 


Chimpanzee 


0.010 


0.003 


3.6* 




Gorilla 


0.009 


0.009 


1.1 




Orangutan 


0.018 


0.020 


0.9 


Chimpanzee 


Gorilla 


0.006 


0.007 


0.8 




Orangutan 


0.014 


0.019 


0.7 


Gorilla 


Orangutan 


0.014 


0.025 


0.6 



The Table 5B ratios were calculated according to Li (1993) J. Mol Evol 36:96; Li 
et ah (1985) Moh Biol. Evol 2:150. For all comparisons, statistical significance was 
calculated by Mests, as suggested in Zhang et al. (1998) Proc. Natl. Acad. Sci. USA 
10 95:3708. Statistically significant comparisons are indicated by one or more asterisks, 

with P values as follows: *, PO.05, **, PO.01, ***, PO.005. Exon sequences are from 
Hacia et ah (1998) Nature Genetics 18:155. GenBank accession numbers: human, 
NM_000058.1, chimpanzee, AF019075, gorilla, AF019076, orangutan, AF019077, 
rhesus, AF019078. 

15 The elevated K A /K S ratios revealed by pairwise comparisons of the human and 

chimpanzee BRCA1 sequences demonstrate the action of positive selection, but such 
comparisons alone do not reveal which of the two genes compared, the human or the 
chimpanzee, has been positively selected. However, if the primate BRCA1 sequences are 
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considered in a proper phylogenetic framework, only those pairwise comparisons which 
include the human gene show ratios greater than one, indicating that only the human gene 
has been positively selected (Table 5B). To confirm that positive selection operated on 
exon 1 1 of BRCA1 exclusively within the human lineage, the statistical test of positive 
5 selection proposed by Zhang et al (1998) Proc. Natl Acad, ScL USA 95:3708-3713, was 
used. This test is especially appropriate when the number of nucleotides is large, as in the 
present case (3423 bp). This procedure first determines nonsynonymous nucleotide 
substitutions per nonsynonymous site (6 N ) and synonymous substitutions per synonymous 
site (bs) for each individual branch of a phylogenetic tree (Zhang et al (1998), supra), 

10 Positive selection is supported only on those branches for which &n - bs can be shown to 
be statistically significant (Zhang et al (1998), supra). For BRCA1, this is true for only 
one branch of the primate tree shown in Figure 9: the branch which leads from the 
human/chimpanzee common ancestor to modern humans, where b^lb% = 3.6. Thus, we 
believe that in the case of the BRCA1 gene, positive selection operated directly and 

15 exclusively on the human lineage. 

While it is formally possible that elevated K A /Ks ratios might reflect some locus 
or chromosomal-specific anomaly (such as suppression of Ks due, for example, to 
isochoric differences in GC content), rather than the effects of positive selection, this is 
unlikely in the present case, for several reasons. First, the estimated Ks values for the 

20 hominoid BRCA1 genes, including human, were compared to those previously estimated 
for other well-studied hominoid loci, including lysozyme (Messier et al (1997), supra) 
and ECP (Zhang et al (1998), supra). There is no evidence for a statistically significant 
difference in these values. This argues against some unusual suppression of Ks in human 
BRCAL Second, examination of GC content (Sueoka, N. in Evolving Genes and Proteins 

25 (eds. Bryson, V. & Vogel, H.J.) 479-496 (Academic Press, NY, 1964)) and codon usage 
patterns (Sharp et al (1988) Nucl Acids Res. 16:8207-821 1) of the primate BRCA1 genes 
shows no significant differences from average mammalian values. 

This demonstration of strong positive selection on the human BRCA1 gene 
constitutes the first molecular support for a theory long advanced by anthropologists. 
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Human infants require, and receive, prolonged periods of post-birth care — longer than in 
any of our close primate relatives. Short, R.V. (1976) Proc. R. Soc. Lond. B 195:3-24, 
first postulated that human females can only furnish such extended care to human infants 
in the context of a long term pair bond with a male partner who provides assistance. The 
5 maintenance of long term pair bonds was strengthened by development of exaggerated 
(as compared to our close primate relatives) human secondary sex characteristics 
including enlarged female breasts (Short (1976), supra). Thus, strong selective pressures 
resulted in development of enlarged human breasts which develop prior to first pregnancy 
and lactation, contrary to the pattern seen in our hominoid relatives (Dixson, A.F. in 

10 Primate Sexuality: Comparative Studies of the Prosimians, Monkeys, Apes and Human 
Beings. 214 (Oxford Univ. Press, Oxford, 1998)). 

Evidence suggests that in addition to its function as a tumor suppressor (Xu et al 
(1999) Mol Ce//3(3):389-395; Shene/a/. (1998) Oncogene 17(24):31 15-3124; Dennis, 
C. (1999) Nature Genetics 22:10; and Xu et al (1999) Nature Genetics 22\11-AZ\ the 

15 BRCA1 protein plays an important role in normal development of breast tissue (Dennis, 
C. (1999), supra; Xu et al. (1999) Nature Genetics 22:37-43; and Thompson et al (1999) 
Nature Genetics 9:444-450), particularly attainment of typical mammary gland and duct 
size (Dennis, C. (1999), supra; and Xu et al (1999) Nature Genetics 22:37-43). These 
facts suggest that positive selection on this gene in humans promoted expansion of the 

20 female human breast, and ultimately, helped promote long term care of dependent human 
infants. This long term dependency of human infants was essential for the development 
and transmission of complex human culture. Because positive selection seems to have 
been concentrated upon exon 1 1 of BRCA1, the prediction follows that the region of the 
BRCA1 protein encoded by exon 1 1 specifically plays a role in normal breast 

25 development. The data provided here suggests that strong selective pressures during 
human evolution led to amino acid replacements in BRCA1 that promoted a unique 
pattern of breast development in human females, which facilitated the evolution of some 
human behaviors. 
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EXAMPLE 15: Characterization of BRCA1 Polynucleotide and Polypeptide 

Having identified evolutionarily significant nucleotide changes in the BRCA1 
gene and corresponding amino acid changes in the BRCA1 protein, the next step is to test 
these molecules in a suitable model system to analyze the functional effect of the 
5 nucleotide and amino acid changes on the model. For example, the human BRCA1 

polynucleotide can be transfected into a cultured host cell such as adipocytes to determine 
its effect on cell growth or replication. 

EXAMPLE 16: Identification of Positively-Selected CD59 

10 Comparative evolutionary analysis of the CD59 genes of several primate species 

has revealed that the chimpanzee CD59 gene has been subjected to positive selection. 
CD59 protein is also known as protectin, lF-5Ag, H19, HRF20, MACIF, MIRL, and P- 
18. CD59 is expressed on all peripheral blood leukocytes and erythrocytes (Meri et al. 
(1996) Biochem. J. 316:923-935). Its function is to restrict lysis of human cells by 

15 complement (Meri et al (1996), supra). More specifically, CD59 acts as one of the 
inhibitors of membrane attack complexes (MACs). MACs are complexes of 20 some 
complement proteins that make hole-like lesions in cell membranes (Meri et al (1996), 
supra). These MACs, in the absence of proper restrictive elements {i.e., CD59 and a few 
other proteins) would destroy host cells as well as invading pathogens. Essentially then, 

20 CD59 protects the cells of the body from the complement arm of its own defense systems 
(Meri et al (1996), supra). The chimpanzee homo log of this protein was examined 
because the human homolog has been implicated in progression to AIDS in infected 
individuals. It has been shown that CD59 is one of the host cell derived proteins that is 
selectively taken up by HIV virions (Frank et al (1996) AIDS 10:1611-1620). 

25 Additionally, it has been shown (Saifuddin et al (1995) J. Exp. Med. 182:501-509) that 
HIV virions which have incorporated host cell CD59 are protected from the action of 
complement. Thus it appears that in humans, HIV uses CD59 to protect itself from attack 
by the victim's immune system, and thus to further the course of infection. 

To obtain primate CD59 cDNA sequences, total RNA was prepared (using either 
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the RNeasy® kit (Qiagen), or the RNAse-free Rapid Total RNA kit (5 Prime - 3 Prime, 
Inc.)) from primate tissues (whole fresh blood from chimpanzees, gorillas, and 
orangutans). mRNA was isolated from total RNA using the Mini-Oligo(dT) Cellulose 
Spin Columns (5 Prime - 3 Prime, Inc.). cDNA was synthesized from mRNA with oligo 
5 dT and/or random priming using the Superscript Preamplification System for First Strand 
cDNA Synthesis (Gibco BRL). The protein-coding region of the primate CD59 gene 
was amplified from cDNA using primers (concentration =100 nmole/|j,l) designed from 
the published human sequence. PCR conditions for CD59 amplification were 94 °C 
initial pre-melt (4 min), followed by 35 cycles of 94 °C (15 sec), 58 °C (1 min 15 sec), 

10 72°C (1 min 15 sec), and a final 72°C extension for 10 minutes. PCR was accomplished 
on a Perkin-Elmer GeneAmp® PCR System 9700 thermocycler, using Ready-to-Go PCR 
beads (Amersham Pharmacia Biotech) in a 50 jil total reaction volume. Appropriately- 
sized products were purified from agarose gels using the QiaQuick Gel Extraction kit 
(Qiagen). Both strands of the amplification products were sequenced directly using the 

15 Big Dye Cycle Sequencing Kit and analyzed on a 373 A DNA sequencer (ABI 
BioSystems). 

As shown in Table 6, all comparisons to the chimpanzee CD59 sequence display 
K A /Ks ratios greater than one, demonstrating that it is the chimpanzee CD59 gene that has 
been positively-selected. 

20 

Table 6: K A /Ks Ratios for Selected Primate CD59 cDNA Sequences 



Genes Compared 


K A /K S Ratios 


Chimpanzee to Human 


1.8 


Chimpanzee to Gorilla 


1.5 


Chimpanzee to Orangutan 


2.3 


Chimpanzee to Green Monkey 


3.0 



EXAMPLE 17: Characterization of CD59 Positively-Selected Sequences 

Proceeding on the hypothesis that strong selection pressure has resulted in 
25 adaptive changes in the chimpanzee CD59 molecule such that disease progression is 
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retarded because the virus is unable to usurp CD59's protective role for itself, it then 
follows that comparisons of the CD59 gene of other closely-related non-human primates 
to the human gene should display K A /Ks ratios less than one for those species that have 
not been confronted by the HIV-1 virus over evolutionary periods. Conversely, all 
5 comparisons to the chimpanzee gene should display K A /Ks ratios greater than one. These 
two tests, taken together, will definitively establish whether the chimpanzee or human 
gene was positively selected. Although the gorilla {Gorilla gorilla) is the closest relative 
to humans and chimpanzees, its postulated historical range in Africa suggests that gorillas 
could have been at some time exposed to the HIV-1 virus. We thus examined the CD59 

10 gene from both the gorilla and the orangutan (Pongo pygmaeus). The latter species, 

confined to Southeast Asia, is unlikely to have been exposed to HIV over an evolutionary 
time frame. The nucleotide sequences of the human and orangutan genes were 
determined by direct sequencing of cDNAs prepared from RNA previously isolated from 
whole fresh blood taken from these two species. 

15 The next step is to determine how chimpanzee CDS 9 contributes to chimpanzee 

resistance to progression to full-blown AIDS using assays of HIV replication in cell 
culture. Human white blood cell lines, transfected with, and expressing, the chimpanzee 
CD59 protein, should display reduced rates of viral replication (using standard assays 
familiar to practitioners of the art) as compared to control lines of untransfected human 

20 cells. In contrast, chimpanzee white blood cell lines expressing human CD59 should 

display increased viral loads as compared to control, untransfected chimpanzee cell lines. 

EXAMPLE 18: Molecular Modeling of CD59 

Modeling of the inferred chimpanzee protein sequence of CD59 upon the known 
25 three-dimensional structure of human (Meri et al 1996 Biochem J. 3 16:923-935) has 
provided additional evidence for the role of this protein in explaining chimpanzee 
resistance to AIDS progression. It has been shown that in human CD59, residue Asn 77 
is the link for the GPI anchor (Meri et al (1996) Biochem 1 316:923-935), which is 
essential for function of the protein. The GPI anchor is responsible for anchoring the 
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protein to the cell membrane (Meri et al (1996), supra). Our sequencing of the 
chimpanzee CD59 gene reveals that the inferred protein structure of chimpanzee CD59 
contains a duplication of the section of the protein that contains the GPI link, i.e., 
NEQLENGG (see Table 7 and Figure 10). 

5 

Table 7: Comparison of Human and Chimpanzee CD59 Amino Acid Sequence 

Human SLQCYNCPNP TADCKTAVNC SSDFDACLIT KAGLQVYNKC 

Chimpanzee SLQCYNCPNP TADCKTAVNC SSDFDACLIT KAGLQVYNKC 

10 Human WKFEHCNFMD VTTRLRENEL TYYCCKKDLC NFNEQLENGG 

Chimpanzee WKLEHCNFKD LTTRLRENEL TYYCCKKDLC NFNEQLENGG 

Human TSLS EKTVLLLVTP FLAAAAWSLHP 

Chimpanzee NEQLENGGNE QLENGG T SLS EKTVLLLVTP FLAAAAWSLHP 



15 



20 



Human (SEQIDNO:12) 
Chimpanzee(SEQ ID NO: 13) 



Italics/underline indicates variation in amino acids. 



This suggests that while the basic function of CD59 is most likely conserved 
between chimpanzee and human, some changes have probably occurred in the orientation 
of the protein with respect to the cell membrane. This may render the chimpanzee 
protein unusable to the HIV virion when it is incorporated by the virion. Alternatively, 

25 the chimpanzee protein may not be subject to incorporation by the HIV virion, in contrast 
to the human CDS 9. Either of these (testable) alternatives would likely mean that in the 
chimpanzee, HIV virions are subject to attack by MAC complexes. This would thus 
reduce amounts of virus available to replicate, and thus contribute to chimpanzee 
resistance to progression to full-blown AIDS. Once these alternatives have been tested to 

30 determine which is correct, then the information can be used to design a therapeutic 

intervention for infected humans that mimics the chimpanzee resistance to progression to 
full-blown AIDS. 
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EXAMPLE 19: Identification of Positively-Selected DC-SIGN 

Comparative evolutionary analyses of DC-SIGN genes of human, chimpanzee and 
gorilla have revealed that the chimpanzee DC-SIGN gene has been subjected to positive 

5 selection. Figures 11-13 (SEQ. ID. NOS. 6-8) show the nucleotide sequences of human, 
chimpanzee and gorilla DC-SIGN genes, respectively. Table 8 provides the K A /Ks values 
calculated by pairwise comparison of the human, chimpanzee and gorilla DC-SIGN 
genes. Note that only those comparisons with chimpanzee show K A /K$ values greater 
than one, indicating that the chimpanzee gene has been positively selected. 

10 Table 8: K A /Ks Ratios for Selected Primate DC-SIGN cDNA Sequences 



Genes Compared 


K A /K S Ratios 


Chimpanzee to Human 


1.3 


Human to Gorilla 


0.87 


Chimpanzee to Gorilla 


1.3 



As discussed herein, DC-SIGN is expressed on dendritic cells and is known to 
provide a mechanism for transport of HIV- 1 virus to the lymph nodes. HIV-l binds to 
the extracellular portion of DC-SIGN and infects the undifferentiated T cells in the lymph 

15 nodes via their CD4 proteins. This expansion in infection ultimately leads to compromise 
of the immune system and subsequently to full-blown AIDS. Interestingly, DC-SIGNS's 
major ligand appears to be ICAM-3. As described herein, chimpanzee ICAM-3 shows 
the highest K A /K S ratio of any known AIDS-related protein. It is not yet clear whether 
positive selection on chimpanzee ICAM-3 was a result of compensatory changes that 

20 allow ICAM-3 to retain its ability to bind to DC-SIGN. 

Although the foregoing invention has been described in some detail by way of 
illustration and example for purposes of clarity and understanding, it will be apparent to 
those of ordinary skill in the art that certain changes and modifications can be practiced. 
25 Therefore, the description and examples should not be construed as limiting the scope of 
the invention, which is delineated by the appended claims. 
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